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We investigate the gravitational property of the quantum vacuum by treating its large energy 
density predicted by quantum field theory seriously and assuming that it does gravitate to obey 
the equivalence principle of general relativity. We find that the quantum vacuum would gravitate 
differently from what people previously thought. The consequence of this difference is an accelerating 
universe with a small Hubble expansion rate H « Ae~? VGA Ls instead of the previous prediction 
H = Jf 8rGipr2? /3 œx VGA? > oo which was unbounded, as the high energy cutoff A is taken to 
infinity. In this sense, at least the “old” cosmological constant problem would be resolved. Moreover, 
it gives the observed slow rate of the accelerating expansion as A is taken to be some large value 
of the order of Planck energy or higher. This result suggests that there is no necessity to introduce 
the cosmological constant, which is required to be fine tuned to an accuracy of 107!°°, or other 
forms of dark energy, which are required to have peculiar negative pressure, to explain the observed 


accelerating expansion of the Universe. 


I. INTRODUCTION 


The two pillars that much of modern physics is based 
on are Quantum Mechanics (QM) and General Relativ- 
ity (GR). QM is the most successful scientific theory in 
history, which has never been found to fail in repetitive 
experiments. GR is also a successful theory which has 
so far managed to survive every test [i]. In particular, 
the last major prediction of GR-the gravitational waves, 
has finally been directly detected on Sept 2015 [2]. How- 
ever, these two theories seem to be incompatible at a 
fundamental level (see e.g. [3]). The unification of both 
theories is a big challenge to modern theoretical physi- 
cists. 


While the test of the combination of QM and GR is 
still difficult in lab, our Universe already provides one 
of the biggest confrontations between both theories: the 
Cosmological Constant Problem [4]. Quantum field the- 
ory (QFT) predicts a huge vacuum energy density from 
various sources. Meanwhile, the equivalence principle of 
GR requires that every form of energy gravitates in the 
same way. When combining these concepts together, it 
is widely supposed that the vacuum energy gravitates 
as a cosmological constant. However, the observed ef- 
fective cosmological constant is so small compared with 
the QFT’s prediction that an unknown bare cosmological 
constant (6) has to cancel this huge contribution from the 
vacuum to better than at least 50 to 120 decimal places! 
It is an extremely difficult fine-tune problem that gets 
even worse when the higher loop corrections are included 


5l. 

In 1998, the discovery of the accelerating expansion 
of the Universe [6] [7] has further strengthened the im- 
portance of this problem. Before this, one only needs 
to worry about the “old” cosmological constant problem 
of explaining why the effective cosmological constant is 
not large. Now, one also has to face the challenge of the 
“new” cosmological constant problem of explaining why 
it has such a specific small value from the observation, 


which is the same order of magnitude as the present mass 
density of the Universe (coincidence problem). 


This problem is widely regarded as one of the major 
obstacles to further progress in fundamental physics (for 
example, see Witten 2001 [8]). Its importance has been 
emphasized by various authors from different aspects. 
For example, it has been described as a “veritable crisis” 
(Weinberg 1989, [4| p.1), an “unexplained puzzle” (Kolb 
and Turner 1993 , [9] p.198), “the most striking prob- 
lem in contemporary physics” (Dolgov 1997 [IO] p.1) and 
even “the mother of all physics problems” , “the worst 
prediction ever” (Susskind 2015 chapter two). While 
it might be possible that people working on a particular 
problem tend to emphasize or even exaggerate its impor- 
tance, those authors all agree that this is a problem that 
needs to be solved, although there is little agreement on 
what is the right direction to find the solution [I2]. 


In this paper, we make a proposal for addressing the 
cosmological constant problem. We treat the divergent 
vacuum energy density predicted by QFT seriously with- 
out trying renormalization and assume that it does grav- 
itate to obey the equivalence principle of GR. We notice 
that the magnitude of the vacuum fluctuation itself also 
fluctuates, which leads to a constantly fluctuating and 
extremely inhomogeneous vacuum energy density. As a 
result, the quantum vacuum gravitates differently from a 
cosmological constant. Instead, at each spatial point, the 
spacetime sourced by the vacuum oscillates alternatively 
between expansion and contraction, and the phases of 
the oscillations at neighboring points are different. In 
this manner of vacuum gravitation, although the grav- 
itational effect produced by the vacuum energy is still 
huge at sufficiently small scales (Planck scale), its effect 
at macroscopic scales is largely canceled. Moreover, due 
to the weak parametric resonance of those oscillations, 
the expansion outweighs contraction a little bit during 
each oscillation. This effect accumulates at sufficiently 
large scales (cosmological scale), resulting in an observ- 
able effect—the slow accelerating expansion of the Uni- 


verse. 

Our proposal harkens back to Wheeler’s spacetime 
foam [13] [4] and suggests that it is this foamy structure 
which leads to the cosmological constant we see today. 

The paper is organized as follows: in section [i] we first 
review several key steps in formulating the cosmological 
constant problem; in section we point out that the 
vacuum energy density is not a constant but is constantly 
fluctuating and extremely inhomogeneous; in section [IV] 
we investigate the differences made by the extreme in- 
homogeneity of the quantum vacuum by introducing a 
simple model; in section [V] we give the solutions to this 
model by solving the Einstein field equations and show 
how metric fluctuations leads to the slow accelerating 
expansion of the Universe; in section we explain the 
meaning of our results; in B WTI investigate the 
back reaction effect of the resulting spacetime on the mat- 
ter fields propagating on it; in section [VIII] we generalize 
our results to more general metrics; in section [X] we dis- 
cuss some questions raised and a couple of new concepts 
suggested by the different way of vacuum gravitation. 

The units and metric signature are set to be c = A = 1 
and (—, +, +, +) throughout except otherwise specified. 


I. THE FORMULATION OF THE 
COSMOLOGICAL CONSTANT PROBLEM 


The cosmological constant problem arises when trying 
to combine GR and QFT to investigate the gravitational 
property of the vacuum: 

G uv + Aoguy = 84GT io >» (1) 
where Gyy = Ruy — 5RQuv is the Einstein tensor and the 
parameter Ap is the bare cosmological constant. 

One crucial step in formulating the cosmological con- 
stant problem is assuming that the vacuum energy den- 
sity is equivalent to a cosmological constant. First, it is 
argued that the vacuum is Lorentz invariant and thus ev- 
ery observer would see the same vacuum. In Minkowski 
spacetime, nuv is the only Lorentz invariant (0, 2) tensor 
up to a constant. Thus the vacuum stress-energy tensor 
must be proportional to nuy (see, e.g. [15], [12]) 


Tiv (t, X) = —p nuv. (2) 


This relation is then straightforwardly generalized to 
curved spacetime: 


Tiv (t,x) = —0™ gu (t,x). (3) 


If Tig“ does take the above form (3), the vacuum energy 
density p’°° has to be a constant, which is the require- 
ment of the conservation of the stress-energy tensor 


VET% = 0, (4) 


The effect of a stress-energy tensor of the form is 
equivalent to that of a cosmological constant, as can be 


seen by moving the term 87GTV7"° in to the left-hand 
side 


Gu =F cf Juv = 0, (5) 
where, 
ret = Av + 8Gp*; (6) 


Or equivalently by moving the term Aggy, in to the 
right-hand side 


Gu = -8TG PXE Juv, (7) 
where, 
Ab 
vac = vac . 8 
Peff P ale 8rG ( ) 


So anything that contributes to the energy density of the 
vacuum acts like a cosmological constant and thus con- 
tributes to the effective cosmological constant. Or equiv- 
alently we can say that the bare cosmological constant 
acts like a source of vacuum energy and thus contributes 
to the total effective vacuum energy density. This equiva- 
lence is the origin of the identification of the cosmological 
constant with the vacuum energy density. 

Following the above formulations, the effective cosmo- 
logical constant Aeg or the total effective vacuum energy 
density p% are the quantities that can be constrained 
and measured through experiments. While solar sys- 
tem and galactic observations have placed a small upper 
bound on Aef, large scale cosmological observations pro- 
vide the most accurate measurement. It is interpreted 
as a form of “dark energy”, which drives the observed 
accelerating expansion of the Universe [6] [7]. 

Based on the assumption of homogeneity and isotropy 
of the Universe, the metric has the cosmology’s standard 
FLRW form, which is, for the spatially flat case, 


ds? = —dt* + a? (t) (dx? + dy? + dz?) . (9) 


Then by applying the equations or (7) for the above 
special metric , one obtains the contributions to the 
Hubble expansion rate H = a@/a and the acceleration of 


the scale factor ä from Agg and/or pYf° are 


3H? = Neg = 8TG PXS, (10) 
Ne 8 G vac 
a= = Avac = 1 ett a. (11) 


The solution to the dynamic equation is 
a(t) = a(0)e™, (12) 


where H is determined by the initial value constraint 
equation (10). 

According to the Lambda-CDM model of the big bang 
cosmology, the effective cosmological constant is respon- 
sible for the accelerating expansion of the Universe as 


shown in and contributes about 69% to the current 
Hubble expansion rate : 


det = 39 Hê ~ 4.32 x 10-**(GeV)?, (13) 
or 
px’ = Oy perit © 2.57 x 10747 (GeV)4, (14) 


where Q, = 0.69 is the dark energy density parameter, 
Ho is the current observed Hubble constant and perit = 
3H 
aC 


is the critical density. 


Unfortunately the predicted energy density of the vac- 
uum from QFT is much larger than this. It receives 
contributions from various sources, including the zero 
point energies (~ 107?(GeV)*) of all fundamental quan- 
tum fields due to vacuum fluctuations, the phase tran- 
sitions due to the spontaneous symmetry breaking of 
electroweak theory (~ 10°(GeV)*) and any other known 
and unknown phase transitions in the early Universe (e.g. 
from chiral symmetry breaking in QCD (~ 10~?(GeV)*), 
grand unification (~ 10°4(GeV)*) etc) [12] [I7]. Each con- 
tribution is larger than the observed value by 50 
to 120 orders of magnitude. There is no mechanism in 
the standard model which suggests any relations between 
the individual contributions, so it is customary to assume 
that the total vacuum energy density is at least as large 
as any of the individual contributions [I2]. One thus has 
to fine tune the unknown bare cosmological constant A» 
to a precision of at least 50 decimal places to cancel the 
excess vacuum energy density. 


1 d kdk’ 1 k-k’ 
e k k’ + 
a | Gar (vi e+ Te 
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because of the terms of the form akap and al al. 


Direct calculation shows the magnitude of the fluctu- 
ation of the vacuum energy density diverges as the same 
order as the energy density itself, 


where 


(21) 


II. THE FLUCTUATING QUANTUM VACUUM 
ENERGY DENSITY 


The vacuum energy density is treated as a constant in 
the usual formulation of the cosmological constant prob- 
lem. While this is true for the expectation value, it is not 
true for the actual energy density. 

That’s because the vacuum is not an eigenstate of the 
local energy density operator Too, although it is an eigen- 
state of the global Hamiltonian operator H = f dèx Too. 
This implies that the total vacuum energy all over the 
space is constant but its density fluctuates at individual 
points. 

To see this more clearly, consider a quantized real 
massless scalar field ¢ in Minkowski spacetime as an ex- 
ample: 


dk 1 —i(wt-k-x) | „t „+i(wt-k:x) 
163) = gega (me e, 
(15) 


where the temporal frequency w and the spatial frequency 


k in are related to each other by w = |k]. 
The vacuum state |0), which is defined as 


ax|0) =0, for all k, (16) 
is an eigenstate of the Hamiltonian operator 


1 
H= f eT = 5 | ero (axa, + ajax) : (17) 
where Too is defined as 
es. 2 
Too = 5 (4? + (V9)?) (18) 


But, |0) is not an eigenstate of the energy density op- 
erator 


C (19) 


? 


where A is the effective QFT’s high energy cutoff. (For 
more details on this calculation, see equation in Ap- 
pendix[A]) Thus the energy density fluctuates as violently 
as its own magnitude. With such huge fluctuations, the 
vacuum energy density p’°° is not a constant in space or 
time. 


Furthermore, the energy density of the vacuum is not 
only not a constant in time at a fixed spatial point, it 
also varies from place to place. In other words, the en- 
ergy density of vacuum is varying wildly at every spatial 
point and the variation is not in phase for different spatial 
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FIG. 1. Plot of the expectation value of the square of the 
energy density difference as a function of spacial separation 
AAg. 


points. This results in an extremely inhomogeneous vac- 
uum. The extreme inhomogeneity can be illustrated by 
directly calculating the expectation value of the square 
of difference between energy density at different spatial 
points, 


({ (Too (t,x) — Too (t, x’) )’\) 


Ap? (Az) = 
Ceo Hot 


» (22) 


where Ax = |x — x’| and we have normalized Ap? by 
dividing its asymptotic value (Too)? (the curly bracket 
{} is the symmetrization operator which is defined by 
(A2)). The behavior of Ap? for the scalar field in 
Minkowski vacuum is plotted in FIG.{1| which shows that 


the magnitude of the energy density difference between 


two spacial points quickly goes up to the order of (Too) 
itself as their distance increases by only the order of 1/A. 
(For more details on the calculations and how the energy 
density fluctuates all over the spacetime, see Appendix 
[A}) 

As the vacuum is clearly not homogeneous, equa- 
tion is not valid as it depends on a homogeneous 
and isotropic matter field and metric. Therefore a new 
method of relating vacuum energy density to the ob- 
served Hubble expansion rate is required. 


IV. DIFFERENCES MADE BY THE 
INHOMOGENEOUS VACUUM—A SIMPLE 
MODEL 


The extreme inhomogeneity of the vacuum means its 
gravitational effect cannot be treated perturbatively, so 
another method is required. As solutions to the fully 
general Einstein equations are difficult to obtain, we will 
first look at a highly simplified model. 

A. Beyond the FLRW metric 


To describe the gravitational property of the inhomo- 
geneous quantum vacuum, we must allow inhomogeneity 
in the metric. This is accomplished by allowing the scale 
factor a(t) in the FLRW metric (9) to have spatial de- 
pendence, 

ds? = —dt? + a° (t, x) (dz? + dy” + dz”). (23) 

The full Einstein field equations for the coordinate 

are 


2 2 2 
a 1 /Va Va 
anaE +3 E -A C) =m an 
2 2 , 2 2 
Giu = —2aä — å? (=) mae 2 (2) Oe eagles (25) 
a a 
Goi = 25 a = re = 81GToi, (26) 
Gy = 29298 BÖN _ seen, 4,9=1,2,8, i#j, (27) 


where V = (01,02, 03) is the ordinary gradient operator 
with respect to the spatial coordinates x,y, z. 


By choosing the above simplest inhomogeneous met- 
ric (23), we are assuming a mini-superspace type model, 
and will choose which of these equations do apply later. 
This treatment might result in inconsistencies as general 
vacuum fluctuations of the matter fields posses rich struc- 
tures that they may not produce spacetime described by 
the metric (23). To fully describe the resulting inho- 


mogeneous spacetime, one needs a more general metric. 
However, as a first approximation, using is relatively 
easy to calculate and leads to interesting results. We are 
also going to do the calculations for a more general metric 


in section [VIII 


B. The fluctuating spacetime 


The role played by the value of vacuum energy density 
in the above equations (24), (25), and is differ- 
ent from (10). The value of vacuum energy density is no 
longer directly related to the Hubble rate H through the 
equation (10). This is evident from the 00 component of 
the Einstein equation (24). The equation is only the 
special case of when the spatial derivatives Va and 
V7a are zero, which requires that the matter distribution 
is strictly homogeneous and isotropic. However, as shown 
in the last section, the quantum vacuum is extremely in- 
homogeneous and necessarily anisotropic, which requires 
Va and V2a be huge. This can be seen through the ij 
component of the Einstein equation (27). In fact, due to 
symmetry properties of the quantum vacuum, we have 
the expectation value of shear stress T;; on the right side 


of (27) 


(Tj) =0, i, j=1,2,3, i479. (28) 

Meanwhile, T;; must fluctuate since the quantum vac- 
uum is not its eigenstate either, and the magnitude of 
the fluctuation is on the same order of the vacuum en- 


ergy density 
(Th) ~ (Too)? - (29) 


This means that the Tj; is constantly fluctuating around 
zero with a huge magnitude of the order of vacuum en- 
ergy density. As a result, in (27), the spatial derivatives 
of a(t,x) must also constantly fluctuate with huge mag- 
nitudes. 

More importantly, since the scale factor a(t, x) is spa- 
tially dependent, the physical distance L between two 
spatial points with comoving coordinates x; and x2 is no 
longer related to their comoving distance Ax = |x; — x9| 
by the simple equation L(t) = a(t)Aa and the observed 
global Hubble rate H is no longer equal to the local Hub- 
ble rate a/a. Instead, the physical distance and the global 
Hubble rate are defined as 


L(t) = [ ya? (t, x)dl (30) 


and 


L [E *tx)VeExadl 


H(t)=2 = (31) 


i JE J@ (t,x) 
where the line element dl = \/dx? + dy? + dz?. 


Equation shows the key difference between the 
gravitational behavior of quantum vacuum predicted by 
the homogeneous FLRW metric (9p and by the inhomo- 
geneous metric (23). 

For the homogeneous metric (9), the scale factor a is 
spatially independent and just reduces to 


H(t) = EA). (32) 


In this case, there are only two distinct choices for Hubble 
rates on a spatial slice t = Const under the initial value 
constraint equation 


a 87Gpvac 
es ee (33) 


which implies that all points in space have to be simul- 
taneously expanding or contracting at the same constant 
rate (Here we do not include the cosmological constant 
A). 
But for the inhomogeneous metric (23), the scale factor 
a is spatially dependent and there is much more freedom 
in choosing different local Hubble rates at different spa- 
tial points of the slice t = Const under the corresponding 
initial value constraint equation (24). 

In fact, the local Hubble rates must be constantly 
changing over spatial directions within very small length 
scales. This can be seen from the initial value constraint 
equations (26), which can be rewritten as 


V (=) = —4nGJ, (34) 


where J = (71, To2, Tog) is vacuum energy flux" 


The solution to or is 


“ (t,x) = Č(t, xo) — 47G J I(t x’) dl, (35) 


where dl’ = (dz’, dy', dz’) and xo is an arbitrary spatial 
point. The above solution shows that the difference 
in the local Hubble rates a/a between xp and x, is de- 
termined by the spatial accumulations (integral) of the 
vacuum energy flux J. Similar to the shear stress, J has 
zero expectation value 


(J) =0 (36) 
but huge fluctuations 
J = (I?) ~ (Too) ~ Af +00, (37) 


which implies that the local Hubble rates differ from 
point to point due to the fluctuations. The average of 
the absolute value of å/a can be estimated with the con- 
straint equation 


f (9 ~ yG (Too) ~ VGA. (38) 


1 One might notice that requires V x J = 0, which means 
that to produce the metric of the form 23), the energy flux of 
the matter field needs to be curl free. As mentioned in the last 
paragraph of section [IV A] this is not true for general matter 
fields, but here as a first approximation we will use to esti- 
mate the magnitude of change in å/a. 


By using (35), we find that the difference in local Hubble 
rates becomes comparable with itself for points separated 


. 1 P 
by only a distance of the order Ax ~ Jar 8 A > +00: 


a (2) ~ 4rGJAr ~ VGN? ~ (Ə); (39) 


Up to this point, we have used the equations (24), 
and (27). These equations are all initial value constraint 
equations which do not contain the scale factor’s second 
order time derivative a. To get the information about 
the time evolution of a(t,x), we also need to use (25). A 
linear combination of equations and gives, 


1 
Goo + — (Gi + G22 + G33) = ——, (40) 
a a 
where all the spatial derivatives of a cancel and only the 
second order time derivative left. Therefore we reach the 
following dynamic evolution equation for a(t, x): 


ä + 07 (t,x)e = 0, (41) 
where 


2 AnG : 1 
Q EEJ p+ > oP; , p= Too, Pi = vg Tu (42) 
i=1 


If Q? > 0, which is true if the matter fields satisfy 
normal energy conditions, describes a harmonic os- 
cillator with time dependent frequency. The most basic 
behavior of a harmonic oscillator is that it oscillates back 
and forth around its equilibrium point, which implies 
that the local Hubble rates å/a are periodically changing 
signs over time. By using equation you can find that 
a/a must also have this periodic sign change in a given 
spatial direction. 

Physically, these fluctuating features of å/a imply that, 
at any instant of time, if the space is expanding in a small 
region, it has to be contracting in neighboring regions; 
and at any spatial point, if the space is expanding now, 
it has to be contracting later. 

These features result in huge cancellations when cal- 
culating the averaged H through (31). The observable 
overall net Hubble rate can be small although the abso- 
lute value of the local Hubble rate |a/a| at each individual 
point has to be huge to satisfy the constraint equation 
(24). In other words, while the instantaneous rates of 
expansion or contraction at a fixed spatial point can be 
large, their effects can be canceled in a way that the phys- 
ical distance would not grow 10!?° times larger than 
what is observed. 

This picture of fluctuating spacetime is not completely 
new. It is similar to the concept of spacetime foam de- 
vised by John Wheeler that in a quantum theory 
of gravity spacetime would have a foamy, jittery nature 
and would consist of many small, ever-changing, regions 
in which spacetime are not definite, but fluctuates. His 


reason for this “foamy” picture is the same as ours—at 
sufficiently small scales the energy of vacuum fluctua- 
tions would be large enough to cause significant depar- 
tures from the smooth spacetime we see at macroscopic 
scales. 

The solution for a(t, x) will be given by equations (51), 
and in the next section[V]to describe this foamy 
structure more precisely. 


C. Methods and assumptions in solving the system 


In principle, we need a full quantum theory of gravity 
to solve the evolution details of this quantum gravita- 
tional system. Unfortunately, no satisfactory theory of 
quantum gravity exists yet. 

In this paper, we are not trying to quantize gravity. 
Instead, we are still keeping the spacetime metric a(t, x) 
as classical, but quantizing the fields propagating on it. 
The key difference from the usual semiclassical gravity 
is that we go one more step—instead of assuming the 
semiclassical Einstein equation, where the curvature of 
the spacetime is sourced by the expectation value of the 
quantum field stress energy tensor, we also take the huge 
fluctuations of the stress energy tensor into account. In 
our method, the sources of gravity are stochastic classical 
fields whose stochastic properties are determined by their 
quantum fluctuations. 

The evolution details of the scale factor a(t,x) de- 
scribed by equation (41) depends on the property of the 
time dependent frequency Q(t, x) given by (42). For both 
simplicity and clarity, in the following sections we inves- 
tigate the properties of Q by considering the contribution 
from a real massless scalar field ¢. In this case, the stress 
energy tensor for a general spacetime metric guv is 


1 
Tw = ViuoVive T 59u V OV. (43) 


Direct calculation using the inhomogeneous metric (23) 
gives that 


3 


i=l 

where all the spatial derivatives and explicit dependence 

on the metric a are canceled. Thus we obtain 

8rG¢? 
3 


= > 0, (45) 
which is not explicitly dependent on the metric a(t, x). 

However, the resulting spacetime sourced by this mass- 
less scalar field ¢ does have back reaction effect on ¢ it- 
self. This is because ¢ obeys the equation of motion in 
curved spacetime 


VV, = (V—99""0,d) =0, (46) 


1 
— Oo 
v/-g * 


which reduces to 
ô, (a°O.¢) — V - (aVe) =0 (47) 


for the special metric (23). 

Incorporating the back reaction effect by solving both 
the Einstein equations (24), (25), (26), for the met- 
ric a and the equation of motion (47) for the field ¢ at the 
same time is difficult. Fortunately, solving the system in 
this way is unnecessary. Physically, the quantum vac- 
uum locally behaves as a huge energy reservoir, so that 
the back reaction effect on it should be small and can 
be neglected. In our method, we will first assume that 
the quantized field ¢ is still taking the flat spacetime 
form of for field modes below the effective QFT’s 
high frequency cutoff A. We use to calculate the 
stochastic property of the time dependent frequency Q 
and then solve to get the resulting curved spacetime 
described by the metric a(t, x). This will be done in the 
next section [V] 

We then investigate the back reaction effect in section 
VIlþy quantizing the field ¢ in the resulting curved space- 
time. It turns out that, while the resulting spacetime is 
fluctuating, the fluctuation happens at scales which are 
much smaller than the length scale 1/A. Therefore the 
corrections to the field modes with frequencies below the 
cutoff A is quite small and thus the flat spacetime quan- 
tization (15) is valid to high precision. (See equations 
(126) (or (158)), (127) and for quantitatively how 
high this precision is.) In this way we justify neglecting 
the aforementioned back reaction. 

Empirically, this must be true since ordinary QFT 
has achieved great successes by assuming flat Minkowski 
background and using the expansion (15). So if our 
method is correct, (15) has to be still valid even the back- 
ground spacetime is no longer flat but wildly fluctuating 
at small scales. In other words, the resulting spacetime 
should still looks like Minkowskian for low frequency field 
modes. Long wavelength fields ride over the Wheeler’s 
foam as if it is not there. This is similar to the behavior 
of very long wavelength water waves which do not notice 
the rapidly fluctuating atomic soup over which they ride. 


V. THE SOLUTION FOR a(t,x) 


In this section we give the solution for the local scale 
factor a(t, x). 


A. Parametric resonance 


One important feature of a harmonic oscillator with 
time dependent frequency is that it may exhibit para- 
metric resonance behavior. 

If the Q(t, x) is strictly periodic in time with a period 
T, the property of the solutions of has been thor- 
oughly studied by Floquet theory [18]. Under certain 


conditions (for example, the condition (78)), the para- 
metric resonance phenomenan occurs and the general so- 
lution of (41) is (see e.g. Eq(27.6) in Chapter V of [19]) 


a(t, x) = cet P, (t, x) + coe x? Py(t, x), (48) 


where Hx > 0, cı and c2 are constants. The P) and P> 
are purely periodic functions of time with period T. They 
are in general functions oscillating around zero. The am- 
plitude of the first term in increases exponentially 
with time while the second term decreases exponentially. 
Therefore the first term will become dominant and the 
solution will approach a pure exponential evolution 


a(t, x) ~ e™t P(t, x), (49) 


where we have absorbed the constant cı into P(t,x) by 
letting P(t, x) = ci Pi (t,x). 

Physically, the exponential evolution of the amplitude 
of a(t,x) is easy to understand. If Q is strictly peri- 
odic, the system will finally reach a steady pattern of 
evolution (when the second term in has been highly 
suppressed). In this pattern, after each period of evo- 
lution of the system, a increases by a fixed ratio, i.e. 
a(t + T,x) = uxa(t,x), which results in the exponential 
increase since after n cycles, a(t + nT,x) = pRa(t,x). 
Here the ux is related to the H, by Hx = Bie, 

Due to the stochastic nature of quantum fluctuations, 
the Q(t, x) in is not strictly periodic. However, its 
behavior is still similar to a periodic function. In fact, 
Q exhibits quasiperiodic behavior in the sense that it 
is always varying around its mean value back and forth 
on an approximately fixed time scale. To see this, we 
calculate the following normalized covariance: 


x (At) = Cov (97(t1, x), 2? (t2, x) (50) 
_ ({(O? (tr) = (07 (t1))) (t2) — (0? (t2))) F) 
((9? — (922))?) 


where At = tı — t2 and we have dropped the label x in 
the second line of the above definition since the final 
result is independent with x. 

Explicit expression for y as a function of At is given by 
(A12), which is plotted in FIG. |2| It describes how 9? at 
different times change around their mean values together. 
We say that two Q? separated by time difference At are 
positively (negatively) correlated if x(At) > 0(< 0), since 
it means that they are most likely to be at the same 
(opposite) side of their mean value (7). 

FIG. [2] and show that Q? at different times are 
strongly correlated at close range. Especially, the nega- 
tive correlation is strongest when At ~ 2/A, which im- 
plies that if at t = 0 the Q? is above its mean value (7), 
then at t ~ 2/A, it is most likely below (Q?). So basi- 
cally, Q? varies around its mean value quasiperiodically 
on the time scale T ~ 1/A. 

This quasiperiodic behavior of Q should also lead to 
parametric resonance behavior seen in (49), instead with 


> 


AAt 


FIG. 2. Plot of the normalized covariance x as a function of 
temporal separation AAt. 


a difference in that H, would become time dependent, 
i.e. the solution would take the following form 


a(t, x) = elo ER P x), (51) 


where P(t, x) here is no longer a strictly periodic func- 
tion as in but a quasiperiodic function with the same 
quasiperiod of the order 1/A as the time dependent fre- 
quency Q(t, x). (The solution (58) for P(t, x) in the next 
subsection reveals this property.) 

The physical mechanism is similar. The system will 
also reach a final steady evolution pattern. In this pat- 
tern, after each quasiperiod of evolution of the system, 
a will increase by an approximately fixed ratio. Suppose 
that during the ith cycle of quasiperiod T;, a increases 
by a factor Hix, i.e. a(t + ft) = pjxa(t,x). Then after 


“so 


Because the quasiperiods T; d the factors Hix e gener- 
ally different from each other, the exponent in (51) would 
need to take the form of integral ion: 

The detailed oscillating behavior of P(t, x) is not ob- 
servable at macroscopic scales. However, the factor of the 
exponential increase e/o #«(t)4t’ can be observed. In fact, 
inserting into (30), the observable physical distance 
would become 


the n cycles, we have a > T;, x) 


L(t) = L(0)e"*, (52) 


where 
=| | \/P2(t,x)dl (53) 


and the global Hubble expansion rate H is 


1 i / / 
= =a H,.(t!)dt!. (54) 


In the next two subsections, we are going to give the 
solution for P(t, x) and the global Hubble expansion rate 
H. 


B. The solution for P(t,x) 


The magnitude of the time dependent frequency 2 is 
of the order ~ \/G (Too) ~ VGA?, while Q itself varies 
roughly with a characteristic frequency A (this has been 
shown by FIG. 2). Then according to (41), the scale 
factor a would oscillate with a period that roughly goes 
as T = 2n/Q ~ 1/VGA? « 1/A, as A > oo, where 
1/A is the time scale on which the Q itself would change 
significantly. 

So comparing to the oscillating period T of the scale 
factor a, the variation of 2 itself is very slow, although 
the time 1/A is already very short for large A. There- 
fore, during one period of the oscillation of a, Q is almost 
constant since it has not have a chance to change signif- 
icantly during such a short time scale. In this sense the 
time dependent frequency 2 is slowly varying and the 
evolution of the scale factor a is an adiabatic process. 

The slow variation of Q can be verified in a more formal 
way by calculating the expectation values of Q? = SG 2 


and (42)” = SzG g2, Using mae we have 


1 
= 377. — * an4 
1 E dk = = GA, (55) 
dQN\?\ 8G 1 si 
a = k 3 
(Op 3 ayi 2” 
_ 8rG 1 5 6 
ee idk = 7 SaN (56) 


(55) just gives Q ~ VGA? as expected, (56) gives 
dQ/dt ~ VGA. Therefore, during one period of oscilla- 
tion T ~ 2r/Q ~ Jan? we have, as A — +00, the slow 
varying condition (see equation (49.1) in Chapter VII of 


20) 
TdQ/dt K Q, (57) 


is satisfied. Thus the system varies adiabatically since Q 
varies only slightly during the one period of oscillation 
time T. 

The leading order solution of the equation (41) for a 
harmonic oscillator with the slowly varying frequency Q 
can be obtained by a first order WKB approximation. 
This adiabatic approximation neglects the small expo- 
nential factor in (51). It gives the solution P(t, x) which 
is describing the oscillating behavior of a(t,x). The re- 
sult is, 


P(t,x) = a cos (f owa + x) . (58) 


The P(t,x) above is a quasiperiodic function with the 
same quasiperiod of the order 1/A as the time dependent 


frequency Q(t, x) just as expected. The two constants of 
integration Ap and 6, in can be determined by the 
initial values a(0,x) and a(0, x). 

The quantum vacuum is fluctuating everywhere, but 
its statistical property must be still the same everywhere. 
Correspondingly, the statistical property of P(t,x) must 
also be the same everywhere, which requires that the 
constant Ao to be independent with respect to the spatial 
coordinate x. In addition, the constant Ag can be chosen 
as any nonzero value since the scale factor a multiplying 
by any nonzero constant describes physically equivalent 
spacetimes. 

The initial phase 6, at different places must be de- 
pendent on x. In applying the initial value constraint 
equation (35), neglecting the small exponential factor in 
(51) and neglecting the relatively small time derivative 
terms of the slowly varying frequency Q, we obtain the 


result, 
4nrG üi 
———— J '). dl’ 
(0, x) f (0, x’) ? 
(59) 


where 6x, is the initial phase of the scale factor a at an 
arbitrary spatial point xo. 

In solutions and we see the fluctuating na- 
ture of spacetime at very small scales as described in 
the previous section [IV B} In particular, shows that 
the phases of a(t,x) vary on a given initial Cauchy slice; 
some locations contract while others expand. In this new 
physical picture the catastrophic vacuum energy density 
is confined to very small scales. 


Q(0, Xo) 


tan Ox, = Q(0, x), 


tan Ox. + 


C. The global Hubble expansion rate H 


As the system is adiabatic, the parametric resonance 
effect is weak. The adiabatic solution in the last 
subsection does not include the parametric resonance and 
thus misses the small exponential factor expected in (51). 
In this subsection we go beyond the adiabatic approxi- 
mation and investigate the exact strength of the weak 
parametric resonance. 

When considering the weak parametric resonance ef- 
fect, the constant Ao in (58) would become time and 
space dependent and take the following form 


A(t,x) = Agelo Hx (tat! (60) 


in order to satisfy (51). 

To determine how the H(t) depends on the spacetime 
dependent frequency Q(t,x), we consider the adiabatic 
invariant of a harmonic oscillator with time dependent 
frequency, which is defined as 


I(t,x) = —, (61) 


where 


E= La + a?) /Q. (62) 


Replace the constant Ao in (58) by A(t,x) and then 
plug it into the above expression (61) we get that 


I(t,x) = TA, (63) 


where we have neglected the time derivatives of A and 
Q in the above equation (63), which are higher order in- 
finitesimals. J is invariant in the first order adiabatic ap- 
proximation. When going to higher orders, J will slowly 
change with time. Through the relation between I 
and A we can obtain how the A(t, x) changes by investi- 
gating how accurately the adiabatic invariant is preserved 
and how it changes with time. 

It has been proved by Robnik and Romanovski [21] [22] 
that, in full generality (no restrictions on the function 
Q(t,x)), the final value of the adiabatic invariant for the 
average energy I = E/Q is always greater or equal to 
the initial value Ig = Eo/Qo (see the references [21] [22] 
for precise definition about the average energy). In 
other words, the average value of the adiabatic invari- 
ant I = E/Q for the mean value of the energy never 
decreases, which is a kind of irreversibility statement. It 
is conserved only for infinitely slow process, i.e. an ideal 
adiabatic process. 

Therefore, in the case of our quasiperiodic frequency 
Q(t,x) in (41), I will also always increase. Moreover, 
it will increase by a fixed factor after each quasiperiod 
of evolution, which results in an exponentially increasing 
I. This is in fact evident because of the weak paramet- 
ric resonance effect. In the following we investigate this 
exponential behavior in detail. 

First we construct the evolution equation for the adi- 
abatic invariant J. Do the canonical transformation 


a = /21/Qsiny, (64) 
å = V2IQ cos ». (65) 


Then the evolution equations for a and its conjugate mo- 
mentum å transfer to the evolution equation for the new 
action variable J and the angle variable y, 


dI Q 

Be Z I” cos? 

Fi g 082%; (66) 
dp Q A 

a~ Q + 59 I 2%. (67) 


Integrating yields 
t 
m= (2 J H(t) , (68) 
0 
where 
Ae) = = 32 (69) 
x(t) = — 56 608 2p. 
The H,(t’) in the above equation (69) is just the same 
with the H(t’) defined in (51) and (60), which can be 
seen by applying equation (63). Thus equation con- 


structed the dependence of H,,(t’) on the time dependent 
frequency Q(t’, x). 


The observable global Hubble expansion rate H is the 
average of H,,(t’) over time, which was defined by equa- 


tion (54). Plugging into gives, 


17-0. 
H = Re (-; arar) : (70) 


When the slow varying condition holds, from equa- 
tion we know that dy/dt is positive, i.e. y is a 
monotonic function in time. Thus we can change the in- 
tegral in from the integration over t’ to integration 
over yy’: 


where yo = y(0) and y = y(t). 

To evaluate H, we formally treat y as a complex vari- 
able and close the contour integral in the upper half 
plane. The integrand in has no singularities for real 
y if the slow varying condition (57) holds. Equation 
implies that y ~ Nt ~ VGA?t, so the length of the in- 
terval y — yo ~ VGA?t goes to infinity as A > +oo. 
Hence the principle contribution to the integral in 
comes from the residue values at singularities y(x) inside 
the contour: 


1 , È siy dt 
H = 7 Re (223k (5e i ew) . (72) 


dy’ 


Each term in gives a contribution containing a 
factor exp (-2 Im P(k))« So the dominant contribution in 
comes from the singularities near the real axis, i.e. 
those with the smallest positive imaginary part. To keep 
the calculation simple, we retain only those terms. Since 
Q(t) varies quasiperiodically with a characteristic time 
T ~ 1/A, the number of singularities near the real axis 
would roughly be on the order t/r ~ At. Therefore the 


H in is roughly 
H ~ Aexp(—2Imyq)). (73) 


Let t(x) be the (complex) “instant” corresponding to 
the singularity yx): P) = Plt) ~ Qt). In general, 
|t()| has the same order of magnitude as the character- 
istic time T ~ 1/A of variation of the Q. Remember that 
Q ~ VGA?, thus the order of magnitude of the exponent 


in (73) is 


Img) ~ OT ~ VGA. (74) 
Therefore, inserting into gives 
H = ahe f VOA, (75) 


where @ and 8 are two dimensionless constants which 
depend on the variation details of the time dependent 
frequency Q(t,x). Therefore H becomes exponentially 
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small in the limit of taking A to infinity. This is a man- 
ifestation of the well-established result that the error in 
adiabatic invariant is exponentially small for analytic Q 
[20] 2I]. In fact, the technique we used in deriving 
is very similar to the one used in deriving the error in 
adiabatic invariant in the pages “160 — 161” of PO]. 


D. A more intuitive explanation 


So far we have obtained our key result for the 
global Hubble expansion rate H. To understand the 
mechanism of weak parametric resonance better, we give 
a more intuitive explanation in this subsection. 

Consider the following simplest parametric oscillator: 


+w’ (t)z =0, (76) 
where 
w? (t) = we (1 + h cos yt). (77) 


The behavior of the above harmonic oscillator with time 
dependent frequency has been thoroughly studied (see 
e.g. eq(27.7) in Chapter V of [I9]). The parametric reso- 
nance occurs when the frequency y with which w(t) varies 
is close to any value 2w9/n, i.e. 


2wo 

Tos (78) 
where n is an integer. The strength of the parametric res- 
onance is strongest if y is nearly twice wo, ie. ifn = 1. As 
n increases to infinity, the strength of the parametric res- 
onance decreases to zero. This is easy to understand since 
as n increases, the varying frequency y of w(t) becomes 
slower compared to the oscillator’s natural frequency wo 
and as n — œ, reduces to an ordinary harmonic os- 
cillator with constant frequency which has no parametric 
resonance behavior. 

Now let us go back to Ce a(t,x). The time 
dependent frequency Q(t, x) in (41) is more complicated 
than the w(t) given in our example (77). However, it can 
be written in a similar form: 


2A 
P, 0) = 05 ( +f dy (f (Y) cos yt + g (7) a) 


(79) 
where 


GA‘ 


6r ” ey 


03 = (0) = 
and f(y), g(7) are operator coefficients, whose exact form 
are given by and in Appendix[A] The behav- 
ior of 0?(t,x) for an arbitrary x is the same with ?(t, 0) 
except phase differences. The power spectrum density of 
the varying part of Q?(t,0) (except for the constant 2 
part) given by (A18) is plotted in FIG. 

Unlike the case (77) where the w(t) varies with a single 
frequency y, the Q(t,0) in varies with frequencies 
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FIG. 3. Plot of the power spectrum density of the varying 
part of 9?(t,0) (except for the constant 3 part). 


continuously distributed in the range (0, 2A) with a peak 
around 1.7A (see FIG. B). From we have that, as 
taking the cutoff frequency A to infinity, Qo ~ VGA? > 
2A. Because of the continuity of the spectrum of 2, we 
can always find integers n such that if 


|G 
oy aa 
n> oq A — +00, (81) 


2 € (0,24). (82) 


then 


So Q(t,x) always contains frequencies 20Q9/n that may 
excite resonances. From we see that n — oo as 
taking the cutoff A to infinity. While as n increases, 
the relative magnitude of the resonance frequency 209/n 
decreases comparing to the a(t,x)’s natural frequency 
Qo. Then for reasons similar to the simplest parametric 
oscillator (76), the strength of the parametric resonance 
of (41) would also decrease to zero. This weak parametric 
resonance effect leads to the global Hubble expansion rate 


H-0, as A--+oo. (83) 


E. Numerical verification 


In this subsection, we do a numerical calculation for 
the evolution equation to verify our result (75), 
which describes the dependence of H on cutoff A. 

In this subsection, Planck units will be used, so all 
instances of Newton’s constant are set to unity, G = 1. 

The main idea is to rewrite the time dependent fre- 
quency Q(t) in phase space. (To see more details about 
this numeric method, please check Appendix [B] Here we 
only list the most crucial results.) For a real massless 
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scalar field, we have 


T 3 Ja ]3 Jo! 
Pahta = Ff Soe 


+ pkpw COS Wt cos wt — 2zkpp'w sin wt cos wt. 


zkzwww sin wt sin wt 


(84) 


This is the Weyl transformation of the operator (22(t). 
Here {x,, Pk} are phase space points of a particular field 
mode with momentum k. Approximately, for a particular 
choice of {£k}, {pk}, we can get an classic equation for 
a: 


a({xx}, {px}, t) + Utah, {px}, t)a({xx}, {px}, t) =0 


(85) 
The observed value a,(t) is the average of a({xx}, {Pk}, t) 
over the Wigner pseudo distribution function 


W({zk}, {pk}, t), which is based on the wave func- 
tion of the quantum field: 


ao(t) = / (1 dndn) al{zr}, {Px}, t)W {ar}, {pr}, t). 
k 


(86) 
If the quantum field is in its ground state, we have 


k 


Weah ht =[] te Fe 87) 
k 


which means {£k}, {pk} are all Gaussian variables. 
Based on this observation, our method to simulate this 
equation is as following: i) at first we generate a set of 
random Gaussian numbers for {£k}, {Pk} ; ii) we solve 
the equation for this particular set of numbers; iii) 
then we repeat the process for another set of random 
numbers until a certain amount of repetitions; iv) The 
result a(t) is the average over all samples we have gener- 
ated. We choose the repetition amount to be big enough 
for the results to converge. The result of a single scalar 
field case is illustrated in Fig. 

We can find that the slope of log |ao(t)| ~ t is decreas- 
ing as we increase the cutoff A as we expect. For the 
single field, Q? (t) = sg? repeatedly reaches zero since 
classically @ is continuous and oscillates from positive 
to negative. Around these zero points the slow varying 
condition for Q? is not satisfied. But because the time 
duration of reaching zero is very short, this would not 
cause the adiabatic expanding scheme to breaks down. 
This point reveals in the numerical calculation. 

The real Universe contains many different quantum 
fields. In Fig. we show the result when we in- 
clude two independent massless scalar fields in which 
Q(t) = (di + $3). In this case Q? would not reach 
zero at almost all times since this can happen only when 
both ġı and @2 pass zero, which is unlikely to happen 
frequently. 

In the two field case, we plot the log( H/A) ~ A graph 
to verify the quantitative relation (75). The result is 
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FIG. 4. Numeric result for log |a.(t)| for a single real mass- 
less scalar field. It shows that as A increases, the slope of 
log |ao(t)| decreases. 
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FIG. 5. Numeric result for log |ao(t)| when two scalar fields 
are present and it shows that as A increases, the slope of 
log |ao(t)| decreases. 


illustrated in Fig[6] We can see that for A > 10, the 
result shows decent linearity, which is what we expected 
since the derivation of (75) is only valid for large A. In 
this case, the two constants a = e46 ~ 100 and 6 = 0.12 
according to the numeric calculation. 


VI. MEANING OF OUR RESULTS 


It is interesting to notice that both and give 
the exponential evolution and predict an accelerated ex- 
panding Universe. However, the underlying mechanisms 
are completely different, which leads to opposite results 
on the predicted magnitude of the observable Hubble ex- 
pansion rate H. 

The solution is based on the assumption that 
quantum vacuum energy density is constant all over the 
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FIG. 6. The plot of log( H/A) over A. The fitting result shows 
that a = e*® ~ 100 and 8 = 0.12 in this two-field case. 


spacetime, which is a necessary requirement if one sup- 
pose that vacuum acts as a cosmological constant. This 
assumption leads to a huge Hubble expansion rate 


7 81Gipre 


3 «x VGA? > +00 


H (88) 
as taking the high energy cutoff A to infinity. 

Our proposal is based on the fact that quan- 
tum vacuum energy density is constantly fluctuating and 
extremely inhomogeneous all over the whole spacetime. 
This fact leads to a small Hubble expansion rate given by 
(75) which goes to zero as taking the high energy cutoff 
A to infinity. 

If we can literally take the cutoff A in to infinity, 
then H = 0. In this sense, at least the “old” cosmological 
constant problem would be resolved. 

In principle, this effective theory is valid only up to a 
large but finite cutoff A, which leads to a tiny but nonzero 
H. Since H > 0 as A > +o, there always exists a 
very large cutoff value of A such that H = /Q)Ho ~ 
1.2 x 10-42 GeV to match the observation, where Ho is 
current observed Hubble constant. 

So our result suggests that there is no necessity to in- 
troduce the cosmological constant, which is required to 
be fine tuned to an accuracy of 107120, or other forms of 
dark energy, which are required to have peculiar negative 
pressure, to explain the observed accelerating expansion 
of the Universe. 

The exact value of A cannot be determined since we do 
not know the values of the two dimensionless parameters 
a and 6 in (75). In principle, we need the knowledge of 
all fundamental fields in the Universe to determine a and 
8, this deserves further investigations in the future and 
might provide some hint on elementary particle physics. 

The value of A should be on the order of Planck energy 
or higher. According to the numerical calculation in the 
last subsection, A ~ 1000Ep if we consider contributions 
to Q? from only two scalar fields. If more fundamental 


fields are included, we expect the value of 8 would in- 
crease and thus decrease the value of A needed. This is 
because that it increases the mean value of 2 and, as a 
consequence, reduces the ratio between the variation of 
Q over its mean value (Q) that Q varies slower. A slower 
Q leads to smaller H since the parametric resonance is 
weaker. 


VII. THE BACK REACTION 


In this section, we investigate the back reaction effect 
by quantizing the field ¢ in the resulting curved space- 
time to justify our method of using the quantized field 
expansion in Minkowski spacetime as an approxima- 
tion. 

The standard way to quantize the scalar field ¢ in a 
generic curved spacetime g,,, is by first defining the fol- 
lowing inner product on a spacelike hypersurface X with 
induced metric h;,; and unit normal vector n” (see e.g. 


23): 
E E [ (610,63 — $30,¢.)n"Vidx, (89) 


where h = det hi; and ¢1, $2 are solutions to the equa- 
tion (46). The above inner product is independent of the 
choice of ©. 

One then choose a complete set of mode solutions uk 
of (46) which are orthonormal in the product (89): 


(uk, uw) = ô(k = k’), (90) 
(uk, uk) = —d(k — k’), (91) 
(Ux, Uz) = 0. (92) 


Ap cos” Od V7¢ 
Q 2 Q 


In order to understand the effect from back reaction, 
we need to find out how the mode solutions of the above 
equation in the resulting curved spacetime change 
from the mode solutions of the equation in the 
flat Minkowski spacetime. 


Physically, the correction to should be small for 
wave modes with frequencies lower than the cutoff fre- 
quency A. That is because the wave length of those field 
modes is larger than 27/A, while our spacetime fluctu- 
ates on the length scale 2r/Q ~ 1/(VGA?) « 2r/A. The 
relatively long wave length modes should not be sensitive 
to what is happening on small scales. This is analogous 
to the situation of sound waves traveling in the medium 
such as air or water or solids. The medium is constantly 
fluctuating at atomic scales, but this fluctuation does not 
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Then the field ¢ may be expanded as 


ọ= 5 (axux + aku) ; (93) 
k 


For the flat Minkowski spacetime, i.e. guv = Nv, (46) 
reduces to the usual wave equation 


$- V? =0. (94) 
In this case, the mode solutions are usually chosen as 


1 1 ; 
—i(wt—k-x) 
(27)3/2 Vu ? e 


where w = |k|. Plugging into just gives the 
usual quantum field expansion (15). 

For our specific metric (23), (46) reduces to (47). In 
this case, since the rate of accelerating expansion is ex- 
tremely small, the back reaction effect due to the macro- 
scopic expansion of the Universe is only important on 
large cosmological time scales. For this reason, we only 
worry about the back reaction due to the wildly fluctu- 
ating spacetime at small scales. i.e. we neglect the small 
exponential factor in and use the form of the a based 
on the solution (58): 


U_(t, x) = 


= „Ai cos x 
a(t, x) ~~ Q(t, x) OS (O(t, )), (96) 
where 
O(t,x) = [ote x0 +60 (97) 
0 


Then (47) becomes 


A/Q i Q 
el ( z co? © + sin 20 $+ (= + tanovo) -Vo=0. (98) 


2 


affect the propagation of the sound wave whose wave- 
length is much larger than the atomic scale. Similarly, 
the propagation of the field modes in the “medium” —the 
spacetime, which is constantly fluctuating on scales much 
smaller than the wavelength of the field modes, should 
also not be affected. 

Mathematical demonstration will be given in the fol- 
lowing subsections. 


A. A simplified toy model 


It is complicated to obtain the mode solutions of 
for a generic stochastic function O(t, x) whose stochastic 
property is determined by the quantum nature of the 


field ¢. To illustrate the underlying physical mechanism 
more clearly, we start with a simplified toy model by 
restricting the phase angle O(t, x) defined by to take 
the following form: 

O(t,x) =0t+K-x, (99) 
where both Q and K are constants and they have the 
same order of magnitude Q ~ |K| ~ /GA?. 

Of course this toy model does not describe the real 
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spacetime sourced by the quantum vacuum since the 2 is 
by no means a constant but always varying, although the 
varying is slow compared to it own magnitude. However, 
this toy model possesses the key property needed — the 
spacetime is constantly fluctuating. It will be convenient 
for visualizing the back reaction effect from a fluctuating 
spacetime. 

After setting the Q = Constant and the phase angle 
O(t,x) to be the form of (99), the equation of motion 


for ¢ becomes 


(1 +cos2 (t+ K-x))$—V7¢— 3Nsin2 (Nt +K-x)¢+tan(Nt+K-x)K-Vo=0, (100) 


2 . 
where we have set Ag = V22 such that the average of the coefficient 40 cos? © before ¢ is 1 for convenience. 


In the flat spacetime case (94), each mode solution u, in contains only one single frequency. However, for 
the above fluctuating spacetime case (100), high frequencies mixes with low frequencies and each mode solution must 
contain multiple frequencies. In fact, since describes a strictly periodic system with time period 7/Q and spatial 
period 7/|K]|, each mode solution ux must change from to the following form: 


+00 
ux (t, x) = eite) Co + 5 cme? ™(QtHK:x) 


m=— oo 


mÆ0 


(101) 


where cm are constants. 


Inserting (101) into (100) and using the orthogonality of e?’”(°+K) | we obtain the following infinite system of 
linear equations: 


m—2 
mth equation: 5 (—1)"*"K - (k + 2nK) cn 
1 2 3 
+ 5 O 2(m—1)Q) 522 (w 2(m—1)Q)—K- (k+2(m—1)K)] cm-1 
+ [(w — 2mQ)? — (k + 2mK)?| cm (102) 
1 
E @- 9st 04 Saw 2(m+1)2)+K-(k+2(m4 1K) ae 
+00 
+ $O (-1y7tr*tk - (k + 2nK) cp 
n=m+2 
=0, m = 0,+1,+2,+3,.. 
In the above calculations, we have used the Fourier series expansion 
-+00 
tanz = —2 Yep sin 2ng (103) 
n=1 


to expand the term tan(Qt + K - x) in (100). 

For the equations of m < —1, we successively add the (m + 1)th equation to the mth equation by the order from 
m = —oco to m = —1; and for the equations of m > 1, we successively add the (m — 1)th equation to the mth equation 
by the order from m = +00 to m = 1. Most terms can be eliminated by these elementary row operations and the 
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above infinite system of linear equations (102) becomes 


ifm < -1, (w — 2(m— 1) Q) (w — (2m + 1) Q) cm-1 


1 
2 
-l5 (w — 2m) (w — (2m +1) M) — (ke + 2mK) - (e+ (2m +1) K) is 


: (w — 2(m + 1)Q) (w — (2m +1) QO) — (k + 2(m+4+ 1)K) - (k+ (2m + 1) K) Cm+1 


$5 w- 2(m+2) N) (w — (2m + 1) 9) cm2 = 0; 


if m= 0, XO (-1)"K- (k+ 2nK) cn 


+(w? — k?) Co 
+l; (w — 20) (+0) +K: (k-+2K)| ĉi 


+00 
+) (-1)"*1K - (k + 2nK) cn = 0; 
n=2 


ifm > 1, 5 (w — 2(m — 2) 9) (w — (2m ~ 1) 2) cm2 


+ K w — 2(m — 1)Q) (w — (2m — 1) Q) — (k + 2(m — 1)K) - (k + (2m — 1) K) Cm—1 
+ K (w — 2mQ) (w — (2m — 1) Q) — (k + 2mK) - (k + (2m — 1) K) Cm 
+ (w — 2(m+1)Q) (w -— (2m — 1) Q) cm+1 = 0. (104) 


To characterize the property of the solutions of this system more clearly, we define the following parameters for 
convenience: 


(105) 


As mentioned before that our effective theory has a cutoff A such that only modes with w,|k| < A are relevant, 


which are much smaller than Q ~ |K| ~ VGA? as A grows large. Therefore, we are only interested in the solutions 
of (102) or (104) when w, |k| <Q, i.e. when €, v — 0. 
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Dividing both sides of (104) by 0? and doing some necessary algebraic manipulations, (104) can be rewritten as 
ifm < —1, 
€ 
ee 
[om —1) a 


3¢ +00 Z n 7 +00 € n 
= 2 ns 2 | + ee 
+| (3 — 26°) m 5 2md Se a) am Am +1) cosy I> (ger) Cm 
3e Tx € ie vU ix € i 
— 2 — — — 2 -H 
+| (3 — 267) (m +1) z 7 2(m+1)6 Lia ami “am Desr +0) (zapr) Cm41 
€ 
+[m+2)- 5] Cm+2 = 0; 
ifm = 0, 
= € e 
5 (—1)” (2nô? + ĝu cosy) cn 4 | 1 +28 + ; Svs + S| c1 + (ê -— v?) co (106) 
€ e Ta 
+|-1 +26? — a ôv cos y + | cı + So(-1)""1 (2nd? + dv cosy) cn = 0; 
n=2 
ifm> 1, 
€ 
[m= 2) = 5] em-2 
3¢ +00 ë n y +00 € n 
— 96? 2 
+| (3 — 267) (m — 1) 5 7 2(m—1)8 ey) any (an 3)6e07 +0) (5,9) Jom 
3¢ +00 € n T +00 A n 
+ | (3 — 267) m 5 2mô 2 in g m Dcosyto) Y (S55) Jen 


+[(m +1)- =| Cm+1 = 0. 


As €,v — 0, the leading order asymptotic solution for {cn} of the above system of linear equations (106) depends 
only on the leading order of the coefficients before {cn}. By keeping only the leading term for each coefficient, we 
obtain that the leading order solution of (106) for {cn} satisfies the following infinite system of linear equations: 


3(3 — 262) -2(3-262) -1 0 


0 0 0 
3 2(3 — 267) —(3 — 267) E 0 0 0 — 0 
0 2 (3 — 267) —3£ — ðv cosy 1 0 0 = 0 
68? —46? —1 +28? e-v? —1 +28? —46? 68? co} = | 0 
0 0 -1 —% — §ucosy 3-26? 2 0 c 0 
0 0 0 = 3—28? 2(3 — 267) 3 vee C2 0 
0 0 0 0 1 2(3 — 267) 3(3 — 262) --- C3 0 
(107) 
I 
We will denote the infinite matrix in the above equa- sion: 
tion (107) by B and its elements by bmn with —oo < To 
m,n < +00. In order to have a nonzero solution, the det(B) = boo Moo + y (=1)” bon Mon, (108) 
determinant of B must be zero. This gives us the disper- A 
sion relation that € and v must satisfy in the asymptotic n#0 
regime €, v — 0. where Mon is the 0,n minor of B, i.e. the infinite de- 


terminant that results from deleting the Oth row and the 
The determinant can be calculated by Laplace expan- nth column of B. Due to the symmetry property of B, 
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we have that, for each n Æ 0, First, we rewrite (107) as the following form: 
bon _ bo,-n, Mon a —Mo,-n: (109) 5 bmnEn = —bmoco; m = 0, 1, 2, 3, no (113) 


0 
which implies that all the terms inside the summation We 


symbol 5> of (108) exactly cancel. Therefore, only the Notice that the matrix elements of B has the following 
108 


first term in (108) survive and thus we have that symmetry properties: 
bmn = —bem—n, ifm,n4~0 (114) 
B) = Mo (8°) (° — v?) = 11 
det( ) 00 (6 ) (e v ) 0, ( 0) bmo = bmo, bon = bo,—n. (115) 
which leads to The above symmetry properties leads to the following 
relation 
e = v?, (111) Cn = —C_n, n # 0, (116) 


which implies that we only need to solve c, for n > 0 to 
solve the whole system. 
For convenience, we define the following new variables 


or equivalently 


w =k’. (112) £n by 

This proves that the usual dispersion relation still holds C= eat, e (i7) 
for low frequency field modes. Then using the relation (116), the infinite system of linear 
After setting e? = v?, we start solving the infinite sys- equations (113) simplifies to the following infinite recur- 

tem (107). rence equations: 

5 3 
(4 — 26°) a1 + 2x2 = z t 60087, (118) 
1 

(3 — 26°) x1 + (3 — 267) 2a + 323 = 5 (119) 
(m — 2)%m—2 + (3 — 28°) (m — 1)£m-1 + (3 — 28°) Mim + (M + l)zm41 =0, if m>3, (120) 


where the dependence on € in the equation (107) or (113) Sequences satisfying (122) must take the following form: 
has been eliminated by introducing the new variables 
En, N Æ 0 through (117) and the solution for £n depends Ym = Deos(md +), m>3, (123) 
only on 6. 

In order to find the general formula for the sequence 


{2m}, we define the following new variables: where D and w are two constants and V is determined by 
mss a 


Ym =(m—1)tm-1+Mim, m23. (121) cos = —1 +8, sind =6/2- 62. (124) 
Then the recurrence equations (120) become 


Combining (123) and (121), the general formula for £m 


Ymn—1 + 2(1 — 8 )ym +Ym41 = 0, m>3. (122) can be obtained by iteration 
| 


mon (OK ynn cos(nd + 4) + (— ian] (125) 
_ - CD a Dsec( sin (CPE 4 TE) sin (SE y ' TE) +22), m> 3. 


Replacing the cm in (101) by £m through (117) we obtain that, as € > 0, the mode solution ux(t, x) is asymptotic 


to 


ux (t, x) = ge) 1+e 5 en Cert Kx) 
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+00 
(126) 


m>=—oco 


m0 


where £m is determined by (116), (117), (118), (119), (120) and (125). 
I 


Using the orthogonality of et?™(9t+K:x), the relative 


magnitude of the correction to uk from the usual plane 
wave mode e~*(“+-k*) in Minkowski spacetime can be 
characterized by applying Parseval’s identity: 


2 


og 
[Ault =e] So a, (127) 
"mo 
From the solution (125) we know that as m — oo, 
1 
En ~ —- (128) 
m 


Thus the summation inside the bracket of (127) con- 
verges and the correction 


[Auk(t, x)| ~e>0, ase 0. (129) 


Thus we have demonstrated that the low frequency wave 
modes (w < A) are almost not affected by the fluctuating 
spacetime with much higher frequency (Q ~ VGA?). 


B. General case 


The methods used and results obtained in the last sub- 
section for the particular simplified toy model can 
be generalized to the generic case (98). To start, we 
rewrite to the following form: 


(1+ fi) ¢— Vo — Qo fod + Kofs-Vo=0, (130) 
where 
2 
fi= Ao cos? © — 1, 131) 
Q 
342 / Q , 
fo= oe (à cos? © + sn) /Qo, 132) 
VQ 
f = | — + tan OVO ) /Ko, 133) 
20, 
% = (2), Ko =(|VO)). 134) 


For convenience, we choose the constant Ap such that the 
average of fi 


(filt,x)) = 0. 


Unlike the toy model (100) we used in the last sub- 
section, (130) is not strictly periodic. However, (130) is 


(135) 


quasiperiodic and its quasiperiod is the same as the pe- 
riod of (100). This property is reflected in the Fourier 
transforms fi(w,k), fo(w,k) and f3(w,k) of the func- 
tions fi (t,x), fo(t,x) and f3(t, x) respectively which are 
defined by 


filt,x) = / dwd?k fi(w, k)et**), (136) 
fo(t,x) = J dwd?k fa(w, k) tt, (137) 
f3(t,x) = l dwd?k fs (w, k)ettt**) (138) 


For the function fı(t, x) defined by (131), after setting 
the constant Ao by (135) and considering the slow vary- 
ing property of Q(t,x) and O(t,x) in both temporal and 
spatial directions, its leading order goes as 

filt, x) ~ cos 20, (139) 


which implies that the Fourier transform fı(w,k) would 
have two peaks centered at 


w = £200, 


|k| = 2Ko. (140) 


For the function f2(t,x) defined by (132), the second 
term which includes the factor sin 20 is dominant since 
the first term which includes the factor 2/0? goes as 
~ 1/A — 0 due to the slow varying condition described 


by and (56). Thus, its leading order goes as 
fo(t,x) ~ 3sin 20, (141) 


which implies that the Fourier transform f2(w,k) would 
also have two peaks centered at 


w = £200, 


|k| = 2Ko. (142) 

Similarly, for the function f3(t, x) defined by (133), the 
second term which includes the factor tan © is dominant 
since the absolute value of the first term which includes 
the factor VQ/(QKo) also goes as ~ 1/A —> 0 due to the 
slow varying property of Q in spatial directions. Thus, 
its leading order goes as 


f3(t,x) ~ tano ZS. 


Then using the Fourier series expansion (103) for tan O, 
we know that the Fourier transform f3(w,k) would have 
infinitely many peaks centered at 


(143) 


w=+2nON, |k|=2nKo, n=1,2,3,---. (144) 


(For a rough calculation of the above Fourier transforms, 
see Appendix 


In addition, we have the zero frequency component (see 
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ilar to the system described by the simplified toy model 
(100). The only difference is that the Fourier transforms 
of the coefficients fı, f2, and fs in (130) spread around 


center points given by (140), (142) and (144) while the 


Fourier transforms of the corresponding coefficients in 
(100) are ideal delta functions exactly located at same 


points given by (140), (142) and (144). 
(130) would take the 


in Appendix |C) 


fiw =0.k=0)~0, i=1,2,3. (145) 


Therefore, the mode solution of 


In summary, the system described by (130) is very sim- form similar to (101): 


u(t, x) = eiet) | og + m dus Bk! up(w', keto" ttk" =) (146) 


k’A0 


where ux(w’,k’) is non-negligible only when w’,k’ are taking values around the centers given by (140), (142) and 
(144). 


Inserting (146) into (130) and replacing the coefficients fı (t, x), fo(t,x) and f3(¢, x) in (130) by the equations (136), 
(137) and (138) and then using the orthogonality of ei(¥'t+k"-x) we obtain the following uncountably infinite system 
of linear equations which are similar to (102): 


(w, k’)th equation : E al) = (k+ K’)’| Ux (w, k’) 
+f du" dk" | (w — (w! = w")? fila K”) — ifo (w — (w! =a fa (W, k”) 


—iKo (k + (K — k”)) - fẹ (w”, k”) | ux (w! — w", K — k”) =0, (147) 


where we have defined the notation ux(0, 0) = coô(0, 0) for convenience. 


To characterize the property of the solutions of this system more clearly, we define the following parameters similar 


to (105) for convenience: 


pl. k’ - fs n k” - fs 
cosu = =—— >, Cosh 


€= —, = —, cosy = —__,, cos U = arr KIBI’ = Wiji] 


Ww |k] Ko k i k’ k ° fs 
V=5> ’ ’ 
Qo Qo Qo |k||k’| |k||f3| 


(148) 
Dividing both sides of (147) by Q gives 


N\A k2 Ik 
(w’, k’)th equation : (« — =) — (v + ort oa, cos 7) ux (w’, k’) 
0 


tal w" 2 il w" 
EO a "g" . "g" 
(<- (5-5) Anei e (F- 5) ) K 
5 Siac |k’| ,_ |k" 
—16 | v cos u OA cos u % 


+f dw" d? k” 


cos x") |f3 (w”, k”) | ug (w — w”, k'— k”) =0. (149) 


Similar to the toy model case, as €, v — 0, the leading order solution of (149) for ux(w’,k’) satisfies the following 
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uncountably infinite system of linear equations: 


if (w’,k’) = (0,0): (€ — v?) 5 (0,0) co 


W334 w" E AE fae now 
T dw’ d°k 0, fi (w ,k )—i Qa f2 (w ,k ) 
0 0 


k” 
+16 ( | cos u") |f3 (w", k”) | uk (—w", =k”) = 0, 


Qo 
if (w’,k’) # (0,0): (—iefo(w’,k’) — idv cos ul|f3 (w, k’)|) co 


[@)-G@)]-e 


, "N2 , " 
oe dw" d k” | E = =) fi (w, k”) +i E _ =) fo (w" k”) 


k” Zk’ 


k’ k” 
iĝ ed cos u’ al cos u” } |f3 (w”, k”) || uk (w — w”, k’ — k”) = 0, (150) 
Qo Qo 


The above uncountably infinite system of linear equations ( ) can also be written formally in matrix form 
similar to (107). We use similar notations that denoting the matrix here by B and its elements by biw k’), (wr, Kr) for 
convenience. 

In order to have nonzero solutions, the determinant of the uncountably infinite matrix B has to be zero, which 
gives the dispersion relations that € and v must be satisfied in the asymptotic region €, v —> 0. 

The “determinant” of B can be formally calculated through Laplace expansion similar to (108): 


where we have used the property (145) in obtaining (150) from (149). 
150 


det B = b(o,0),(0,0)M(0,0),(0,0) +f dw" Pk" (—1) eb 0), (ur see) 0,0), (wt Je”) (151) 


w £0 
k” 40 


where M,0,0),(w",k”) is the (0,0), (w”, k”) minor of B, i.e. the ‘determinant’ resulting from deleting the (0,0)th row 
and (w”,k”)th column of B. 
Notice that since fı(t, x), fo(t,x) and f3(¢,x) are all real, their Fourier transforms fı(w,k), fo(w,k) and f3(w,k) 


defined by (136), (137) and (138) must satisfy the following relations: 
fi(w,k) = fi(—w, -k)*, fo(w,k) = Ja(—w, —k)*, f3(w, k) = f;(-w, —k)*, (152) 


where the * means complex conjugate. 
The above symmetry property (152) leads to 


6(0,0), (wr k) = 8(0,0),(—-w",-k)), Moo) (wk) = —M0,0),(-w,-k), i£ (—w”, —k”) # (0,0), (153) 


which implies that all the terms inside the integral symbol f of (151) exactly cancel. Therefore, only the first term 
in (151) survives and thus we have 


det B = Mio,0),(0,0)(€? — v?) = 0, (154) 
which gives again the usual dispersion relation 
e=v o wee. (155) 
After setting the dispersion relation ma we only need to solve the (w’,k’) 4 (0,0)th equations in since 
det B = 0 implies that the (w’,k’) = (0,0)th equation is redundant. 


For convenience, we define new variables x(w’, k’) similar to the z,, defined in (117): 


ux(w’,k’) = ecory(w’,k’), (w’,k’) Æ (0,0). (156) 


Then (150) can be rewritten as 
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if (w/,K’) # (0,0): (B) 2 o ar (w, K’) 


, 
d "k" w L 
+f ut OS | G 


[k"| 


„< ( (k'| ; 
(K COS p 


=i fo(w’,k’) + id cos u|f3 (w, k’)]. 


(157) 


Replacing the ux, (w’,k’) in (146) by z_(w’,k’) through (156) we obtain that, as e + 0, the mode solution ux (t,x) 


is asymptotic to 


u(t, x) = “ges 1+ ae dw! Bk! a (w', kje ttx) l 
w 
k'#0 


where £(w’, k’) is determined by (157). 


Analogous to (125) in the simplified toy model, 
2, (w’, k’) would also go as 


1 
t,(w", k’) ee at (159) 
m 
when w’, k’ taking values around the centers 
w ~+2mO9, |k’|~ 2mKo, m=1,2,3,--- (160) 


(x_~(w’,k’) is negligible if w’,k’ is far away from these 
centers). 

Due to Parseval’s theorem, implies that the inte- 
gral inside the bracket of converges which is similar 
to (127) and thus the correction to u(t, x) also goes as 
€. 

Therefore, when we quantize the scalar field ¢ in our 
wildly fluctuating spacetime by expanding it in terms 
of the annihilation and creation operators according to 
(93), the leading order would still be the form of the 
Minkowski quantum field expansion (15). The correc- 
tion to the dispersion relation w? = k? and the plane 
wave mode e~'(“t-*-*) are on the order ~ e. In addition, 
the extra wave modes which mixing in or 
are all modes with frequencies higher than No ~ VGA?, 
which is much larger than our effective QFT’s cutoff A. 
These extremely high frequency modes beyond the cut- 
off are irrelevant to our low energy physics. This also 
explains why the ordinary QFT works by assuming fixed 
Minkowski spacetime. The small scale structure averages 
out in its effect on the long wavelength low energy fields. 

In summary, we have argued that although our space- 
time sourced by the quantum vacuum is highly curved 
and wildly fluctuating, the back reaction of the resulting 
spacetime on the quantum field sitting on it is small. This 
justifies our method of neglecting back reaction and using 
the quantum field expansion in Minkowski spacetime 
at the beginning. 


(158) 


VIII. THE MORE GENERAL 
INHOMOGENEOUS METRICS 


In previous sections we assume the simplest inhomo- 
geneous metric to describe the spacetime resulting 
from the inhomogeneous vacuum. In this section, we try 
to generalize the result to more general inhomogeneous 
metrics. 

The quantum fluctuations of the vacuum is not com- 
pletely arbitrary, the magnitude of the fluctuations must 
be the same everywhere and in every direction, i.e. 
the spacetime is still stochastically homogeneous and 
isotropic. Thus we can always choose a special gauge and 
construct the following general synchronous coordinate: 

ds? = —dt? + has(t,x)dx"dx’, a,b=1,2,3. (161) 

For the above metric (161), we employ the initial value 
formulation of general relativity. In this formulation, the 
Einstein equation is equivalent to six equations for the 
evolution of the second fundamental form 


kay = — RÈ) — (trk) Kay + 2kack€ 


1 (162) 
+4rGpha, + 8TG (Ta — hate ) j 
plus the usual four constraint equations, 
RC) + (trk)? — kapk™ = 167Gp, (163) 
Dak; = D,(trk) = 8TG jo, (164) 


where kap = than, k® = BPN ha, trk = h®kap, 
p = Too, jo = hłToa, rT = h®Ta, R® is the 3- 
dimensional spatial curvature and D, is the derivative 
operator associated with hap. 


Taking trace on both sides of and then combining 

with (163) gives: 
h? kay — kark = —47G (p + trT). (165) 

It is interesting to notice that there are no spatial deriva- 
tives included on the left hand of the above equation 
(165). The key evolution equation (41) for a(t,x) we 
used in previous sections is just the special case of the 
above equation (165). 

Direct calculation using the expression shows that, 
the contribution from a real massless scalar field to the 
right-hand side of is 


p+trT = 2¢°, (166) 
where all the spatial derivatives of ¢ and all the explicit 
dependence on the metric g,, in the definition of stress 
energy tensor are canceled. It is also interesting to 
notice that the above exact expression is exactly 
the same with the corresponding expression (44) for the 
simplest inhomogeneous metric case. 
We first consider the following special case: 


a° (t,x) 0 0 
halt, x) = 0 tx) 0 (167) 
0 0 @(t,x) 


The spacetime described by the above coordinate 
possesses more freedoms than and thus would ex- 
hibit richer structures. In this case, the expansion rate at 
the same point becomes directionally dependent. Along 
the three principle axes ĉ, y and 2, which are eigenvec- 
tors of the symmetric matrix hab in (167), the expansion 
rates a/a, b/b and é/c can be different. This means that, 
at one same point, the space can be expanding in one or 
two directions and contracting on the other two or one 
directions. 


Under the coordinate system (167), equation (165) be- 
comes 


424E = on 
a b c 


(168) 
which is a generalization of the key evolution equation 
(41) we used in previous sections. 


Let 


ä b 
A = -Q$ (t,x), 37 


ë 
b -93 (t, x), c = -93 (t,x), 
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then (168) immediately leads to 


3 (t,x) + 02 (t,x) + Q(t, x) = 47G (p +trT). (170) 
As the functions a, b and c are alternately symmetric, 
their expectation values must be equal 
(Os) = ae (p+trT), i=1,2,3. (171) 
Unlike equation (41), here Q? (t, x) does not necessarily 
go exactly the same as 42° zG (p +trT). Since 4rG(p+trT) 
is slowly varying, 0? mast also be slowly varying func- 
tions. Otherwise we would have two or three fast varying 
functions sum together and precisely cancel each other to 
give a slowly varying function, which is almost impossible 
in the system with such huge quantum fluctuations. Thus 
the evolution of a, b and c are still adiabatic processes 
and the conclusion we obtained in previous sections still 
holds. We can obtain solutions similar to that 


a(t, x) ~ el m= Pi (t,x), (172) 
b(t, x) = efa Falt) P(t x), (173) 
e(t, x) = efo 3x0 P(t, x), (174) 
and on average 
H; =aAe YCA i=1,2,3, (175) 
where 
1 ; / f 
h= i H;x(t)dt (176) 
0 
Therefore, the observable physical volume is, 
w= [vit t,x) = er. (177) 
where h = det hap = a7b?c?. 
Next we investigate the most general case 
a?(t,x) d(t,x) e(t,x) 
hav(t,x) = | d(t,x) b°(t,x) f(t,x) (178) 
2 


e(t,x) f(t,x) (t,x) 

In this case, the three orthogonal eigenvectors of the 
symmetric matrix hap can rotate in space. This gives 
more freedom and structure to the spacetime evolution 
than in the case described by the coordinate system 
(167). For example, an initial sphere will distort toward 
an ellipsoid with principle axes given by eigenvectors of 
hab, with rates given by time derivatives \;/A; of the 
corresponding eigenvalues \?(t, x), i = 1,2, 3. 

Expanding the dynamic equation using the met- 


(169) ric (178) gives 
| 
@hiyä  Whinbd Chee dhi} d ehisë | fhis f 
h a h b h c h d h e h f 


+ F (hav, hab) = —47G (p + trT), (179) 


where h = det(hap) is the determinant of the matrix 
(178), hž, is the matrix’s (a, b) cofactor and F is a nonlin- 
ear function of the metric components hap and their first 
time derivatives hab. (179) is difficult to handle. Further 
investigations are needed in the future. 

However, the results we obtained for the coordinate 
(167) suggest that, for the most general case (178), the 
eigenvalues A?(t,x) should also evolve adiabatically sim- 
ilar to a?, b? and c°. In other words, we expect that the 
results (172), (173), and can be generalized 
to A; in the most general case and the physical volume 
of space would expand as 


V(t) =| Vh(t, x)d x 


= / Me edad T 
=V (0)e?**, (180) 


where H is determined by (75). 


IX. DISCUSSION 


So far, we have presented a new mechanism of vacuum 
gravitation and showed that it leads to a slow acceler- 
ating expansion instead of a catastrophic explosion of 
the Universe. In this section, we discuss some questions 
raised and a couple of new concepts suggested by this 
different way of vacuum gravitating. 


A. Lorentz invariant cutoffs 


A potential concern is that the cutoff A we are using 
is not Lorentz invariant. However, the results would not 
change if using Lorentz invariant cutoffs instead. This is 
because the A is just used for comparing the magnitude of 
different infinities, whose leading order dependencies on 
A can also be obtained just by dimensional analysis. In 
previous sections, we have taken this simple non-Lorentz 
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invariant cutoff A for convenience. In this subsection, 
we use more complicated but Lorentz invariant Pauli- 
Villars type cutoffs to show directly that the results do 
not change. 

First, we calculate the two point function 


dwd?k i T 
(d(t1, x) b(t2, x) = J (27) w2 — k2 + E ar 
(181) 
where At = tı — t2. We then replace the photon propa- 
gator 


— 182 
w2 — k? + ie (182) 


in (181) by 
1 1 
w? — k? +ie w2—k?— A? pie’ (1853) 


where A^; can be thought of as a fictitious heavy photon, 


which can serve as a Lorentz invariant cutoff. 
Then (181) becomes 


dwdďk _, 1 1 j 
_ N2 a. d 
‘ af (2Q7)4 i 0 “(kR aA? +iej ’ 


(184) 
where we have used the identity 
1 . da 
= ; 185 
AB | (A+(B-A)a)’ a 


(184) is still logarithmically divergent when setting 
At = 0. To make it converge, we employ another Pauli- 
Villars type cutoff Ao by replacing the 


1 
T 
w? — k? — ay, A? + ie ee) 
in (184) with 
1 1 
w? —k?—ayA?+ie w2?—k? —a, A} — AZ ic 
(187) 


Then (184) becomes 


1 


dwdk _; 1 1 : 
_ aas f ka. d T | dal, 
iATAS (Ont e ; ai : daz o &2 (w? — k? — ay A? — ag A2 + ie)? (w? — k? — a A? — a AB + ie)?’ 


(188) 


where we have used the identity (185) again in obtaining (188). 


Performing the integration f d?k first, (188) becomes 


A2A4 1 1 1 , 
1 zJ aa, f dos f dof, f duvet 
167? Jo 0 0 
1 1 


(189) 
1 


l 2_ A2 a 2_ Ar , T ~ = a 
w? — A (ai, aa) +ie yw? — A?(a4, a5) + ie (V -Elora + ie ya Relan ag) tie) 
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where 
A? (a1, a2) = a1 A? + a2 A2, Al? (a, ah) = aA? + ab AZ. (190) 


Then let w = A(a1,a2)u, (189) becomes 


ae dos f dos f dogr | duet a2)Atu T (u, a1, a2, a), (191) 
167? A4 (ay, a2) i 
where the integrand is, 
Io(u, a1, @2,05) (192) 
1 1 1 


= zi 
af — 1 +z = 2 Al? (01,05) | __ie 2 A’? (01,04) i 
Tan az) y" A2(a1,02) | A2(a1,02) je -1+ RT, 4 Klana) i CAET 
If ag 4 a, the integrand is multivalued in complex plane and has four branch points 


Al ie 
U1,2,3,4 = = 1 Ss ae < 1 ie (193) 


Assuming At > 0, the integral contour goes around the lower half i Without loss of ae assuming 


+ 
a2 > ah, we can choose the branch cut being the line connecting u1 = 4/1 — ae € and uz = x Se 
1 


re u — U1,3, we can easily see that the 
U—-U1,3 i 


The integral 
Io is determined by the integral along this branch cut. Since Io > 


integral along the branch cut converges to a finite number. 


ie 
A2 (a1, a2) 


—5/2 
If ag = a4, Io becomes (u -1+ ) /8. The integral can be solved easily in this case, i.e. ho dulo = + 
when e€ > 0 , which is also finite. 
Therefore, after setting At = 0, we have 
_ APAS 
g) T dor f doa f ETE Jim (u,a4,02,05)~ A?, as An Ay ~ A> +o. (194) 
~~ 161? “(ay, Q2) 


Similarly, we have 


; d d 
(#) = Jim, ai gg Otxa) 


AM 
af don f dos f dea, an) fave To(u,a1, 02,05) ~ Af, as A~ Ay ~ Ag +o, (195) 
167 CRE 


which gives the same result as (55). 


When calculating the quantity (3?) from (191) using sharp cutoff A and on the Lorentz invariant Pauli-Villars 


the same technique as in (195), we would get the integral type cutoffs A, Aa, (As, Aa) ae the ee As we men- 
4 neq i : ; í tioned in the beginning of this subsection, these results 
f duu*Io, which is logarithmically divergent. To make it f : 
eee are natural since they could have been guessed by dimen- 
converge, we can employ another two Pauli-Villars type 


: . sional analysis. Therefore, the slow varying condition 
Sint and ha ph ae wit woe Thom (186) to (L87). (57) is still satisfied by using Lorentz invariant cutoffs 


The exact dependence of (8) on Ai, Az, As and A4 iS that our results would not change. 
complicated, but the result for its leading order goes as 


(3?) ~ AS, (196) B. The singularities at a(t,x) =0 


which gives the same result as (56). 


So far, we have obtained that the leading order depen- Four way OF vardum GE Witting, the apace slice 


; 5 P tively expanding and contracting at each spatial point, 
dencies of (ġ ) and (ð ) on the non-Lorentz invariant and, during each such cycle, the expansion outweighs 


the contraction a little bit due to the weak parametric 
resonance effect. This process gives a slowly increasing 
amplitude A(t,x) of the scale factor a(t,x), whose ob- 
servable effect is just the accelerating expansion of our 
Universe. 


Probably one of the biggest concerns about this phys- 
ical picture is the appearance of singularities at points 
a(t, x) = 0—according to the solution (58), the scale fac- 
tor a(t, x) must go through zero whenever the space at x 
switches from contraction phase to expansion phase. 


Singularities are a generic feature of the solution of 
Einstein field equations under rather general energy con- 
ditions (e.g. strong, weak, dominant etc.), which is guar- 
anteed by Penrose-Hawking singularity theorems [24H29]. 
In this paper, since we investigate the gravitational prop- 
erty of quantum vacuum without modifying either QFT 
or GR, the appearance of singularities is inevitable— 
QFT predicts a huge vacuum energy, and according to 
GR, huge energy must collapse to form singularity. 


It is usually thought that the Einstein field equations 
break down at singularities and thus the spacetime evolu- 
tion will stop once the singularity is formed. However, it 
is not the case for our solution to the key dynamic evolu- 
tion equation (41). describes the oscillating motion 
a harmonic oscillator. It is natural for a harmonic os- 
cillator to pass its equilibrium point a(t,x) = 0 at max- 
imum speed without stopping. So in our solution, the 
singularity immediately disappears after it forms and the 
spacetime continues to evolve without stopping. Singu- 
larities just serve as the turning points at which the space 
switches from contraction phase to expansion phase. 


1. Resolving singularity by multiplying a 


In order to understand why in our solution the singu- 
larity is not the end of spacetime evolution, it is helpful to 
review one crucial step in deriving from (40). Rigor- 
ously speaking, we can only obtain the following equation 


from (40): 


a 


- = = (t,x), (197) 


more step—multiply both sides of (197) by a. 


which is not equivalent to (41). To get (41), we need one 


u(t, x) = cge *t-*™) (: + ie 


At the singularities (198), the terms sin 2m(Qt + K - x) 
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Mathematically, a is not allowed to be zero in 
since it is in the denominator. In fact, when writing down 
the Einstein field equations (24), (25), and (27), it 
has been presumed that a 4 0 since if a = 0, the metric 
would become degenerate (g = det(g,,) = —a°® = 0), the 
curvature would become infinite and the Einstein tensor 
are simply not defined there. 

But, after the inequivalent algebraic manipulation of 
multiplying both sides of by a, a is allowed to evolve 
to zero in the resulting equation since there is noth- 
ing wrong for a harmonic oscillator to go through its equi- 
librium point. In this sense, we have smoothly extended 
the solution beyond the singularity by the mathematical 
operation of multiplying both sides of by a (or more 
generally by some power of the metric determinant). 

The idea of resolving a singularity by mulptiplying Ein- 
stein equations with some power of the determinant of 
the metric is not new. Einstein himself had proposed 
this idea with his collaborator Rosen in 1935 (for which 
they credited this idea to Mayer) [80]. Ashtekar used a 
similar trick in his method of “new variables” to develop 
an equivalent Hamiltonian formulation of GR [BI]. It is 
also proposed by Stoica that the equations obtained after 
multiplying the usual Einstein equations by some power 
of the metric determinant are actually more fundamental 
than the usual Einstein equations [32}H40]. In this sense, 
we argue that our spacetime with singularities due to 
the metric becoming degenerate (a = 0) is a legitimate 
solution of GR. 


2. Singularities do not cause problems 


While singularities are natural and inevitable in solu- 
tions to Einstein’s equations, we must discuss the conse- 
quences they bring to this calculation. 

Will the singularities cause serious problems? At least 
in our case we do not feel they cause problems. To see 
this, we investigate how the singularities affect the prop- 
agation of the field modes in our toy model (100). 

In this toy model, the spacetime have singularities at 
the hypersurfaces 


1 
OQt+K-x=(n+ zT n = 0, +1, +2, +3,--- (198) 
Using the relation £m = —£-m, which is evident from 


(116) and (117), the asymptotic mode solution (126) be- 


comes 


+00 
XO tm sin 2m (Nt + K - ») (199) 
m=1 
| 
of (199) are all zero and thus we have 
u(t, x) = cget-**) (200) 


So ux is normal at singularities which shows that the field 
can naturally pass the singularities without problems. 

One might still worry about the divergences of the time 
or spatial derivatives of ux, at the singularities (198). 
However, these divergences arise from those small high 
frequency corrections (terms inside the summation sym- 
bol X of (199)) with frequencies 2mQ which are much 
higher than our effective QFT’s cutoff A. When looking 
at low energy scales (< A), uk behaves the same as the 
mode solution when the background spacetime is 
flat; only when looking at high energy scales (> VGA?) 
which are much higher than the cutoff scale A, those 
small high frequency corrections are noticeable. 

In this sense, the singularities do not cause problems 
at the observable low energy regime—after all, the singu- 
larities only appear (and immediately disappear) above 
Planck energy scales, which should not affect the low en- 
ergy physics whose energy scale is far below Planck. 


C. Similarity of effects of vacuum energy in 
non-gravitational system and gravitational system 


Vacuum fluctuations and their associated vacuum en- 
ergies are direct consequences of the Heisenberg’s uncer- 
tainty principle of quantum mechanics. Although it is 
still controversial [41], various observable effects are of- 
ten ascribed to the existence of vacuum energies and have 
been experimentally verified, which strongly suggests the 
reality of vacuum fluctuations. These vacuum fluctuation 
effects include the spontaneous emission [42], the Lamb 
shift [43], the anomalous magnetic moment of the elec- 
tron and the Casimir effect [46H49]. The reality 
of the vacuum energy associated to the spontaneous sym- 
metry breaking of electroweak theory has also been con- 
firmed by the discovery of the Higgs boson at the LHC 
[50] 51]. 

If we assume that the vacuum fluctuations do exist as 
evidenced by the above listed observable effects, then ac- 
cording to the equivalence principle, the associated vac- 
uum energies would gravitate as well as all other forms 
of energy. This has been experimentally demonstrated 
by, for example, the gravitational test of Lamb shift en- 
ergy [52H54]. The gravitational property of Casimir en- 
ergy has not been tested experimentally, but has been 
demonstrated theoretically with the conclusion that it 
does gravitate according to equivalence principle [55H58]. 

However, in the literature, the value of vacuum energy 
density is usually thought to play a different role in non- 
gravitational systems and in gravitational systems. The 
actual value of the vacuum energy density is generally re- 
garded as irrelevant in non-gravitational contexts based 
on the argument that only energy differences from the 
vacuum are measurable; while when gravity is present, 
the actual value of the energy matters, not just the dif- 
ferences, since the source for the gravitational field is the 
entire energy momentum tensor that its large value may 
be potentially disastrous. 
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We argue differently in this section with the following 
points: (i) the value of vacuum energy density can also be 
relevant in non-gravitational contexts; (ii) the huge value 
of vacuum energy density is not a direct observable and 
that it is not disastrous in a theory of gravity. Moreover, 
there is essentially no difference between the roles played 
by vacuum energy in non-gravitational systems and in 
gravitational systems. In other words, although tech- 
nically more complicated when gravity is included, the 
gravitational effect of the vacuum energy on spacetime 
metric is intrinsically the same as its effect on material 
bodies when gravity is excluded. 


1. Value of vacuum energy is relevant in Casimir effect 


Let us first consider the Casimir effect. The Casimir 
force is usually derived by calculating the change in vac- 
uum energy due to the presence of the conducting plates, 
which acts as mirrors to reflect electromagnetic waves 
(We will call them mirrors in the following). This deriva- 
tion is straightforward, but loses some important physical 
details about what is going on in the system [59}|60]. Due 
to quantum fluctuations, the zero point fields constantly 
impinge on both sides of the mirror and then reflect back, 
which transmit momentum to the mirror and thus re- 
sult in forces on both sides of the mirror. The Casimir 
stress (force per unit area) is just the difference between 
the pressure exerted by the electromagnetic field vacuum 
from inside and outside 


S(t, z, y) = pos = outside, (201) 


where we have set that the two parallel mirrors are nor- 
mal to the z axis. Since the vacuum fluctuations be- 
tween the two mirrors are different from the vacuum fluc- 
tuations outside, the expectation values of TÌ»®side and 
Toutside Would be different and thus gives a net aver- 
age force. Although both Tiside and Teutside are diver- 
gent, this average force is finite since the quartic diver- 
gent Minkowski zero point fluctuations are canceled after 
the subtraction in and one obtains the well known 
Casimir stress 


r2 


(3) = ~aaogt 


Thus the effect of the value of zero point energy disap- 
pears in the calculations. It is for this reason that al- 
though the Casimir effect is usually regarded as evidence 
of the reality of zero point energy, the actual value of its 
energy density is thought to be irrelevant in this effect. 
However, the value of zero point energy density does 
have an effect. Note that only gives the expectation 
value of the Casimir stress S, but S is never a constant, 
it fluctuates. That’s because the amount of momentum 
carried by the zero point fields which impinge on both 
sides of the mirror is constantly fluctuating due to the 
fact that the vacuum is not an eigenstate of the zz com- 
ponent of the stress energy tensor T,,. The magnitude 


(202) 


of the fluctuation of each T,, is large and diverges as the 
same order of the vacuum energy density (Too). For a 
perfect mirror, since the fields on the two sides fluctuate 
independently of each other, the mean-squared stresses 
on the two sides simply add, resulting in the magnitude 
of the fluctuation of the net stress also diverges as 


(AS?) = ((S = (8))?) ~ (Too)? > œ. (203) 


For more realistic imperfect mirrors which become trans- 
parent for frequencies higher than its plasma frequency 
A, the (Too) in contains contributions only from 
field modes of frequencies lower than A and the mean 
squared value of the net stress S goes as 


(AS?) ~ (To)? ~ AB. 


The plasma frequency A in acts as an effective cut- 
off which depends on the microstructure of the mirror. It 
is similar but distinct from the effective QFT’s cutoff A in 
(75), which depends on the microstructure of spacetime. 

Therefore, the value of zero point energy density is 
still physically significant even in non-gravitational sys- 
tem. Its value appears in and (204) to characterize 
the strength of Casimir stress fluctuation, which implies 
that the net Casimir stress is constantly fluctuating with 
huge magnitudes around its small mean value (202). Due 
to this huge fluctuation, at almost any instant, the mag- 
nitude of the stress at each single point of the mirror is 
as large as the value of the zero point energy density. 

However, this effect is strong only at small scales. Its 
measurable effect becomes small at larger scales. In prac- 
tice, the measurements must be taken over some finite 
time interval T and some finite surface area of order /?. 
More precisely, what the force detector measures is the 
time and surface average 


(204) 


S= | dtacdyf(t,x,y)S(t2,y), (205) 
where the averaging function f satisfies 
J dtavayf(t,2, y)=1. (206) 


The exact shape of the averaging function depends on 
the measuring apparatus. On physical grounds one can 
choose f to be a single peak over a time interval T com- 
parable to the experimental resolving time and over a 
spatial region of area 1? comparable to the resolution of 
the measuring device. Although the magnitude of the 
fluctuations of the net stress S is formally infinite as 
shown in (203), the magnitude of the measurable fluc- 
tuations of its average S$ is finite. This is because the 
effect of the vacuum fluctuations at small scales is signif- 
icantly weakened when averaging over larger scales. The 
calculations have been done by G Barton in [61] with the 
conclusion that, for the realistic case where | < cT, the 
mean squared deviation 


(a53) = ((58-(8))")=FE, (207) 
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where the “constant” here is a pure number as could have 
been foreseen on dimensional grounds. The above equa- 
tion shows that (AS?) increases as T decreases, 
which means that the better the measuring device, the 
stronger fluctuation due to the effect of the value of 
the zero point energy density can be measured. And 
in principle, using a perfect instantaneous measuring de- 
vice (T + 0), one can measure the infinite fluctuations of 
the Casimir stress on a perfect mirror due to the infinite 
value of zero point energy density. In practice, however, 
(AS?) is too small to be measured for a real force detec- 
tor whose resolving time T is too large [6]]. 


2. Effect of vacuum energy on the motion of mirrors 


The value of zero point energy density also has effects 
on the dynamic motion of small material bodies. Imag- 
ine that we place a single mirror of very small size in the 
vacuum and then release it. The mirror would experience 
a fluctuating force exerted by the quantum field vacuum 
and starts to move. The equation of motion of the mir- 
ror, which is called quantum Langevin equation, can be 
generally described by 

K =F (t,%,%,6,4,...), (208) 
where X is the mirror’s position, ¢ represents the field 
interacting with the mirror which is usually taken to be 
a scalar field for simplicity and we have set the mirror’s 
mass M = 1 for convenience. The average force in this 
case would be zero because of symmetry 


(FY =0, (209) 
and similar to the Casimir stress fluctuation (203), the 
force here also undergoes wild fluctuations with a mag- 
nitude 
(F?) œ (Too)? > 00. (210) 
The mathematically infinite fluctuating force F gives 
infinite instantaneous accelerations of the mirror through 
(208). Similar to the case of infinite Casimir stress fluc- 
tuation (203), this infinite fluctuating force and infinite 
instantaneous acceleration make sense since they are also 
only significant at very small scales and will not result in 
infinite fluctuation of the mirror’s position at observable 
larger scales. In fact, the mirror would oscillate back and 
forth with very high speeds, but its range of motion is 


still small [62H66]. 


More precisely, suppose that the mirror is initially lo- 
cated at X(0) = 0 with velocity X(0) = 0 and is then 
released at t = 0. The magnitude of its acceleration 


X(t) and velocity X(t), which can be characterized by 
the quantity (X7()) and (aY is large. But, the 
magnitude of the range of the mirror’s fluctuating mo- 


tion, which can be characterized by the observable mean 
squared displacement (X?(t)}, is still small. 


In this sense, the value of vacuum energy density is still 
relevant even in non-gravitational physics. This value 
appears in the equation to characterize the strength 
of the force fluctuations acting on the mirror at small 
scales and it may have small observable effects at larger 
scales such as diffusions predicted in [64H66]. 


3. Analogies between the motion of mirror and the motion 
of a(t, x) 


Although technically more complicated in gravity, the 
basic dynamic equation of motion satisfied by the 
scale factor a(t, x) is in fact very similar to the equation 
of motion satisfied by the mirror’s position X(t). 
Consider only the contribution from the massless scalar 


field ø, (41) is just the following same form as the equa- 
tion (208) 


i = F (a,4), (211) 
where 
F (a, ô) = -FC ga, (212) 


Also, the average of the fluctuating force F (a, 6) is zero 


due to symmetry 


(F (a, $) ) =0, (213) 
and its magnitude of fluctuation 
(F° (a, $) ) x (Tho)? > o0. (214) 


The above two statistical properties and sat- 
isfied by the “force” driving the “motion” of the scale 
factor a are the same with the statistical properties 
and satisfied by the force driving the motion of the 
mirror. In this sense, the role played by the value of the 
vacuum energy density in gravitational system is similar 
to its role in non-gravitational system. 

Concretely speaking, the vacuum energy density re- 
sults in large instantaneous acceleration X and velocity 
X of the mirror, but the observable position fluctuations 
of the mirror, which can be characterized by the quan- 
tity (X o7 is not large. Analogously, the vacuum energy 
density results in the large instantaneous “acceleration” 
ä and “velocity” å of the scale factor, but the observable 
physical distance defined by (30), whose value is deter- 
mined by the quantity (a?), is also not large. These prop- 
erties about a(t,x) are evident from the solutions (51), 
and (75), from which we can see that the quantities 
(a?) and (a?) are as large as (Too)? and (Too) respec- 
tively, while the magnitude of the quantity (a?) is on the 
order 1. 

In this sense, the role played by vacuum energy in 
gravitational system is similar to its role in the non- 
gravitational mirror systems—it appears both at 
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and to show the strongness of vacuum fluctuations 
at microscopic scales (for mirrors, microscopic means 
atomic scale; for gravity, microscopic means Planck scale) 
and their observable effects are both small at macroscopic 
scales. 

By this same kind of mechanism, the violent gravita- 
tional effect produced by the vacuum energy density is 
confined to Planck scales, and its effect at macroscopic 
scales—the accelerating expansion of the Universe, due to 
the weak parametric resonance is so small that, it is only 
observable after accumulations on the largest scale—the 
cosmological scale. 
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Appendix A: Real Massless Scalar Field 


In this appendix we give the calculation details about 
how the quantum vacuum fluctuates all over the space- 
time by using the massless scalar field as an example. 

We first define the covariance of the energy density 
operator at two spacetime points x = (t,x) and a2’ = 
(v’,x’) 


Cov (Too(£), Too(2’)) 
=({ (Too(#) — (Too(x)) ) (Too(x") — (Too(2’)) ) }) (AI) 


where the curly bracket {} in (A1) is the symmetrization 
operator which is defined as, for any two operators A and 
B 


bi 


{AB} = ; (AB + BA). 


Inserting and into gives the following 


result 


(A2) 


Cov (Too(£), Too(2’)) 

1 f BkdPk! wo! +k- k)? 

52) r w2 (A3) 
- COS (Ww +w')At—(k+k’)- Ax), 


where At = t — t and Ax = x — x’ are time and space 
separation of the two spacetime points x and 2’. 


If x and xv’ are timelikely separated, we can find a refer- 
ence frame to set Ax = |Ax| = 0. In this case, evaluation 


1 


Cov (Too(x), Too(2’)) = TETE 


( [—(AAt)® + 21(AAt)* — 72(AAt)? + 36] cos(2AAt) 
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of the integral in (A3) for a high frequency cutoff |k| = A 
gives 


(A4) 


+6 [(AAt)” — 8(AAt)? + 12AAz] sin(2AAt) + 12 [(AAt)? — 6AAt] sin( AAt) + 36 [(AAt)? — 2] cos(AAt) + s) ; 


If x and z’ are spacelikely separated, we can find a reference frame to set At = 0. In this case, evaluation of the 


integral in (A11) for a high frequency cutoff |k| = A gives 


1 


Cov (Too(x), Too(2")) = AF 


+16 [(AAz)? — 6] cos(AAzx) — 64AAz sin(AAz) + 63) 


As At and Az goes to 0, both (A4) and (A5) reduces 
to the variance of the energy density, 


( (oo — (Too) )’) = : (A) = ; (Too)? (A6) 


We then investigate the Pearson product-moment cor- 
relation coefficient 


_ Cov (Too(x), Too(2’)) 


Pax! = } 
OxOg' 


(A7) 


where 


oe = y| Tola) — (Tolo) P). (48) 


The correlation coefficient pz. shows by its magnitude 
the strength of correlation between two random variables. 
Px,x’ is positive if the energy density Too at x and x’ are 
most possibly lying on the same side of the vacuum ex- 
pectation value (Too) = A*/(1677). Thus a positive cor- 
relation coefficient pz, implies the energy density at x 
and x’ tend to be simultaneously greater than, or simul- 
taneously less than the expectation value. Similarly, a 
negative pr, implies the energy density tend to lie on 
opposite sides of the expectation value. We will call the 
energy density Too at x and x’ are positively correlated 
if Px,» > 0 or negatively correlated (anticorrelation) if 
Prax! <0. 

Because of transnational invariance, px, is only de- 
pendent on the temporal and spatial separation At = 
t — t/, Ax = x —x’. For the real massless scalar field 
(15), the behavior of the correlation coefficient py." as 
a function of temporal separation AAt for the case of 
Aa = 0 and as a function of spatial separation AA for 
the case of At = 0 are plotted in FIG. |7| and [8] respec- 
tively. 


( [2(AAz)* — 34(AAz)? + 33] cos(2AAzx) — [12(AAz)? — 50A Az] sin(2AAz) 


(A5) 
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FIG. 7. Plot of correlation coefficient pv s as a function of 
time separation AAt in the case Ax = 0. 


In the temporal direction, i.e. the case of Ax = 0 (Fig. 
F), the correlation coefficient goes quickly from 1 down to 
around —0.9 in a time scale around At = 1.9/A and then 
goes up to 0.7 in a time scale around At = 3.8/A and 
then goes down and up alternatively from positive values 
to negative values with decreasing amplitudes. It roughly 
oscillates as — cos(2AAt)/(AAt)? with a period 7/A as 
At is large. Thus at the extremely small time scales 
At ~ 1.9/A, (A > +00), the energy density are strongly 
anticorrelated. In other words, if at some time the value 
of the energy density is larger than its expectation value, 
for example, by an amount of 0.82 (Too), after a short 
time At = 1.9/A, its value is most likely to be smaller 
than the expectation value, for example, by an amount 
of 0.74 (Too). The difference is 1.56 (Too) only after such 
a short time. 

In the spatial direction, i.e. the case of At = 0 (Fig. 
[8), the correlation coefficient goes quickly from 1 down to 
around —0.14 in a length scale around Ag = 3.24/A and 


AAX 


FIG. 8. Plot of correlation coefficient pz, as a function of 
spatial separation AAz in the case At = 0. 


then goes up to 0.03 in a length scale around Az = 5.4/A 
and then goes down and up alternatively from posi- 
tive values to negative values with decreasing ampli- 
tudes. Compared to the temporal direction, the decay 
in the oscillation amplitude of the correlation coefficient 
is faster in spatial direction. It roughly oscillates as 
2cos(2AAz)/(AAx)* with a period 7/A as Az is large. 
These properties show that the strength of the correla- 
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tion between energy densities at close range in spatial 
direction is not as strong as in the temporal direction. 
For larger spatial separations, p,,./ approaches zero and 
the vacuum energy density Too at different x and a’ fluc- 
tuate independently. These properties result in extreme 
spatial inhomogeneities of the quantum vacuum which 
can be characterized by the quantity Ap? defined by 
in section [II 
The quantity Ap? is related to pz, by 
Ap? = I= prx. (A9) 
The behavior of Ap? has been plotted in FIG. 
Next we calculate the x(At) defined by in section 


Wick expansion of gives 


x(At) = 


2 (s(t, x) am 


where the correlation function can be calculated directly 


by inserting (15) 
: ‘ 1 A ; 
(6(ti.)6(to9) = gly f Be Pak 
An? Jo 


Plugging (A11) into (A10) gives the following result 


(A11) 


x(At) = — z (36 (—2 + AP At?) cos( AAt) + (36 — 72A7At? + 21A4At* — ASAt®) cos(2AAt) 


+ 6 (6 + 2AAt (—6 + A? At?) sin( AAt) + AAt (12 — 8A? Ad? + A*At*) sin(2AAt)) ) À 


The behavior of x(At) has been plotted in FIG. 
of time difference At in the case Ar = 0 (FIG.|7). 
Next we derive the equation in section 


02(t,0) = sL d kdk Vw 
wA  (2r)3 2 


+ (axa — ala.) cos (w+w')t+i (axa — akak) sin (w + w’) r| À 


Specially, the vacuum state |0) is an eigenstate of the 
operator coefficients of the first two terms in the above 
expression (A13p. If k Æ k’, the eigenvalues of the op- 
erator coefficients of the first two terms are zero. Thus 
in this case, the first two terms have to both take zero 
values. If k = k’, the second term is zero since in this 
case w = w’ and thus the factor sin(w—w’)t = 0. So only 
the first term survives and gives the expectation value of 
0?(t, 0): 


a2 = (0 aa 


Ëk w _ GA 
(2r)? 2 6r” 


(A14) 


(A12) 


It is closely related to the correlation coefficient py,’ as a function 


First, 2?(t,0) can be expanded as 


=| (axa, + ajax) cos(w—w')t+i ( akal, | ala) sin (w — w’) t 


(A13) 


For the operator coefficients of the last two terms in the 
expression (A13), the vacuum state |0) is not an eigen- 
state. So the last two terms are constantly fluctuating, 
and the time varying of Q? comes from these two terms. 


After some algebraic manipulations, (A13) can be 
rewritten as the form of for the vacuum state |0), 
where 


161? Bkk Vww! 


ySwtw'<yt+dy (2T)? 2 


( AKA! 


(A15) 


167? 


kdk Vww' . 
g(y)dy = a 


I coven 

(A16) 

Evaluating the above integrals gives the expectation 
values 

(f(y)dy) = (gy)dy) = 0, (A17) 


(Oan?) = (i)a?) = a (40 - 140% + 168 (4) -P(A 


The above expression (A18) gives the power spectrum 
density of the varying part of Q?(t,0) (except the con- 
stant 2 part), which has been plotted in FIG. 


Appendix B: Wigner-Weyl] Description of Quantum 
Mechanics and Numeric simulations 


This chapter explain the principle of the numeric cal- 
culations in the main text. Same as the numeric part 
in the main text, we set G = 1 in this section. Wigner 
functions and Weyl transforms of operators offer a for- 
mulation of quantum mechanics that is equivalent to the 
standard approach given by the Schrödinger equation. 
The Wigner distribution function is a quasi distribution 
function in the phase space. For a particular quantum 
wave function 7(x), its Wigner function is defined as 
(B1) 


Wp) = f dye (e+ Dura -¥) 


The Weyl transform of an quantum operator A is defined 
as 
Amp) = f dye Ele- 4) (B2) 
Then the expectation value of the operator A under the 
state y(x) can be written as 
(4) = | facdpw(e,n)Alc,») (3) 
These two transformations give the Wigner-Weyl discrip- 
tion for quantum mechanics. The expectation values of 
physical quantities are obtained by averaging their Weyl 
transforms over phase space. 
For a harmonic oscillator with frequency w and m = 1, 
the ground state Wigner function is a Gaussian distribu- 
tion function for both x and p 


p? 2 


1 
Wo(x, p) = zE ° (B4) 


We can easily check that the Weyl transform of an op- 
erator H(%) ( or H(p)) is simply replaced the operator 


EE (axa ~ aba 
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i their fluctuations 


Los yo A; 


A18 
ifA<y> 2A. ( ) 


& by x (or p by p). Other than that, another particular 
transform we are going to use in this write-up is 


(7 
2 2 (Ba 
We can see that the transform of the product does not 
necessarily equal to the product of transforms. In the 
following part we are going to get the general expression 
for the transform of the product. 
Before that we notice that Weyl transform can be used 
to construct the original operator , i.e. 


A wena TTY | ép(e—y) 
(lâl) = = f dpa nye 


Using this formula we can construct the transform of 
product of two states: 


ĉp — rp +=; px > xp 


(B6) 


[eve + 5 ABle — e 


= f dzdyle + ¥Alz)(2|Ble— Bye 


dzdydpı dpze'?! (t+ 3—*) e—tP2(®#— 3-2) e—tpy 


~ An? 
z +y/2+z zr—y/2+z 
“Al ul Pi) B( ul , p2) 
= dz dzodp,dpge'7! (P27P) e172 (P-P1) 
Ar? 
Z z 
A(x + >. pi)B(a + >, p2) 
(B7) 
Here we define 
a= Tte (B8) 
a =- +z-r (B9) 


We Taylor-expand A(x + %4, p1) and B(a+%, p2) around 
x and have 


Alz + Fim) = D AM ep) /2)” (B10) 
n=0 ` 
Blo + 3p) = Y BO (w.p2)(e2/2)" (B11) 


and use the facts 


dzz” et?” = (—i)"6™ (y) (B12) 


Qn 
and 


J dy” (y) f(y) = (C10) (B13) 


Therefore, we can write the Weyl transform of operator 
AB as 


Ds antl 


n,m 


CH) ap) Be (ap) (B14) 


The generalized FRW scale factor a satisfies the equation 


4+07(t)a=0 (B15) 
in which 
ea = STga) (B16) 


Now we replace all the quantities by operators, assuming 
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For a massless real scalar field, we can write it as 


: Bk. Pe 
d= 1 amar (ts cos(w,t) + n” sin(wzt)) (B19) 
in which 
1 
Êk = TE + bk) (B20) 
QW 
Pr = TE — by) (B21) 


are the generalized ĉ p operators for each field modes. 


We can write the Weyl transformation of the ((t) 


8 d kd? k! 
O({re}, {pe}, t)? = = | One CELE WWE SIN Wet Sin Wet 
+ PkPk' COS Wet COS Wet — 2LkPk'Wk SIN Wet COS Wy t 
(B22) 


that operators still satisfy the previous equation 
a+ Q(t)?4 =0 (B17) 
This expression is quadratic in x, and px, so if we apply 
with it to (B14), only m+n < 2 terms survive. Assuming 
Bad a({xz}, {pk}, t) is the Weyl transform of operator â, we 
Q? (t) = 3° (t) (B18) have the equation for a as 
| 
E R i 00? ða NQ? ða 1 oN? Oa 0? 0? Ora 070? Oa, 
a+ a H = 
2 ke Ork opk opk Ork 8 kk! OX pKOL pK! Opp ODE: OprODp: OxpOLK Oxr,OpE! OppOX pK! 


The observed value a is the average over Wigner function 


W ({xx}, {pe}, t) 


ao(t) = J (11 nut) a({rx}, {pr}; t)W ({xx}, {Pk}; t) 
k 

(B24) 

If the quantum field is in the ground state, then by 


1 Se 
War} {r} t) = [] <e “ee 


k 


(B25) 


a. Local approximation Generally the equation 
depends on not only the value of Q and a on a par- 
ticular phase space point (x, p), but also on the neighbor- 
ing values (i.e. derivatives). If our solution a is ”smooth” 
enough in the phase space then we can neglect the last 
two derivative terms in the (B23). It can be simplified to 


&+07%a=0. (B26) 


(B23) 


Assuming the length of the Universe is L. We can replace 
the integral by summations. For simplicity, we define 


ta = (B27) 
L 
In > 4 /2W2rn T 27n (B28) 
L L 
D : B29 
Bu 4] 5 Paes (B29) 
L 
The equation can be written as 
4 52 
a+ gpa) a=0 (B30) 


with 


Qi)? = x Von! (Zač z sin nt sin n't 
fin! 
+ Pap, cos nt cos n't 


— £4 P,~ sin nt cosn’t)) (B31) 


= 5 Vn(žasinnt — pa cos nt) 
P 


Here 7 = (n1, n2, n3), n1,2,3 € Z and n = |ñ]. {Xz} {pr} 
are random Gaussian variables with unit standard devia- 
tion. We can solve the equation for a randomly generated 
set of {%;} and {pz}, and repeat. The result a(t) is the 
average over all solutions as long as our sample size is big 
enough. 


Appendix C: Fourier transforms of the coefficients 
in (130) 


In this appendix, we demonstrate the property of the 


spectrum of the coefficients in given by (140), 


and (144). Observing that the cos 20, sin20 and tan O 


in (139), (141) and (143) respectively can all be decom- 


posed as Fourier series sum of the form e’?”°, where 
n = +1,+2,---, we only need to analyze the spectrum 
of ei2n0 


For simplicity, we only analyze the time component 
Fourier transform of e’?"°. The spatial part has similar 
property. The phase angle O is determined by Q through 
while 2 is determined by (79). The power spectrum 
of Q? is given by (illustrated in FIG. |3). 

Calculation of the Fourier transform of e”"© exactly 
based on is complicated. For simplicity, we assume 
that Q taking the following simple form which is similar 


to (77) 


Q = Qo(1 + h cos yt), (C1) 
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where y take the peak value of the power spectrum 
which is around ~ 1.7A (see FIG. |3) and h < 1 to make 
sure that Q > 0. 

Then we have 


hQ 
© = Not + — sin yt. (C2) 
Y 
Using the Jacobi-Anger expansion we have 
= 2nhQ 
ei2n0 Z Do Jml n 0) eino +m )t (C3) 


m=—oo Y 
where Jm is the mth Bessel function of the first kind. 
As |m| — oo, we have 


Anh 1 /nhQ,\'™! 


m! i (e 
which drops faster than the exponential. Therefore, the 
Fourier transform of e’?”© is centered around 2nQo. 

To estimate the magnitude of the Fourier coefficients 
of e'2”© around zero frequency, we evaluate the Bessel 
function for 


m ~ —2n/y ~ VGA > oo. (C5) 
In this case, the zero component Fourier coefficient is 
asymptotic to (see [67]) 


en (v—tanh v)|m| 


4/27|m| tanh v 


since v is determined by h = sechy < 1 that we always 
have v — tanh v > 0. 

When calculating the Fourier transform of e ex- 
actly based on (79), the spectrum becomes continuous 
instead of discrete. But the distribution of the spectrum 
should be similar. 


|Jm(—hm)| ~ > 0, (C6) 
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